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t 

This  report  was  prepared  for  the  Naval  Air  Systems  Command,  Washington, 
D.C.,  under  contracts  N000 19-8 l-C-0395  and  N00019-84-C-0123 ,  "Computer  Code  for 
Flutter-Critical  External-Store  Configurations".  Funding  was  provided  via 
j  Dr.  Daniel  Mulville,  AIR-310B.  The  contract  technical  monitor  was 

Mr.  George  Maggos,  AIR-5302C. 

The  report  consists  of  three  volumes.  Volume  I,  entitled  "User's  Manual", 

I  rovides  instructions  for  using  the  ESP  program  and  presents  descriptions  of 

typical  output.  Volume  II,  "Final  Report  on  Program  Enhancement  and  Delivery", 
Describes  the  work  that  was  performed  under  the  two  contracts.  A  listing  from 
CDC  compilation  of  the  program  is  contained  in  Volume  III,  "Program 
,'ompi  1  at  ion"  . 

The  contributions  of  many  individuals  to  the  successful  completion  of  the 
i  .i tracts  are  gratefully  acknowledged.  Ms.  Ann  Marie  Novak  performed  much  of 
-.lie  work  required  to  convert  the  original  IBM  code  to  a  CDC  version.  Highly 
-i  i  cable  consulting  support  was  provided  by  Mr.  Richard  Chipman,  the  primary 
Developer  of  the  original  ESP  version,  and  by  Mr.  Dino  George  and  Dr.  Joel 
Markowitz,  key  developers  of  FASTOP.  Assistance  on  computing  problems  was 
.  ided  by  several  persons  at  Grumman,  including  (in  alphabetical  order)  Mr. 
•narles  Bores,  Mrs.  Linda  Ehlinger,  Mr.  Joel  Halpert,  Mr.  Luke  Kraner,  Mr. 
or  ild  Mackenzie,  Mr.  Mario  Mistretta,  Mr.  John  Ortgiesen,  Ms.  Florence 
tnpfneimer ,  and  Mrs.  Noreen  Wolt.  Key  contributions  to  making  the  ESP 
/■•a gram  operational  on  the  NADC  Central  Computing  System  were  made  by  Messrs, 
-pert  Richey  and  Howard  Ireland  of  the  Naval  Air  Development  Center. 

' ....  i  1  v ,  Mr.  Louis  Mitchell  of  the  Naval  Air  Systems  Command  provided  valuable 
i:'  ight  icto  program  features  which  would  be  important  to  practicing  flutter 

and  also  provided  helpful  suggestions  during  the  preparation  of  this 

•  rf  . 
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1  -  SUMMARY 


A  pilot  computer  program  for  determining  flutter-critical  external-store 
configurations  has  been  developed  and  made  operational  on  the  Naval  Air 
Development  Center  Central  Computer  System,  The  new  program,  designated  F.SP 
(External-Stores  Program),  is  a  derivative  of  the  previously  developed  Flutter 
And  STrength  Optimization  Program  (FASTOP) , 

Three  key  ingredients  in  ESP  are: 

o  A  p-k  flutter-solution  algorithm  that  includes  an  automated  procedure 
for  defining  the  flutter  speed. 

o  A  calculation  of  the  derivatives  of  the  flutter  speed  with  respect  to 
the  store  parameters. 

o  A  gradient-directed  numerical-search  algorithm  to  seek  out  progres¬ 
sively  lower  flutter  speeds  within  a  parameter  space  defined  by  the 
range  of  store  parameters  at  each  aircraft  store  station. 

A  restart  capability  is  available  for  continuing  a  search  following  a 
visual  review  of  the  search  status  after  a  run.  Also,  an  analysis-only  option 
permits  a  user  to  avoid  the  inclusion  of  input  data  required  only  for  perform- 
*’\s  a  st ore -parameter  search.  To  facilitate  the  introduction  of  data  required 
it  i.e  dynamic  idealization  of  an  aircraft,  an  automated  interface  has  been 
neve  loped  between  ESP  and  both  the  COSMIC  and  the  MacNea 1 -Schwendler 
C'”-porat i on  versions  of  NASTRAN . 

The  version  of  ESP  that  exists  as  of  this  writing  was  developed  primarily 
h  rh-»  obiect  Lve  of  quickly  evaluating  the  feasibility  of  the  store-search 
noepi.  .  ";!ie  program  structure  has  nor  been  optimized  for  minimum  computing 

not  .  11  options  in  the  program  have  been  checked.  Thus,  the  current 
!  i  v r ra I  on  shim  Ld  be  considered  as  a  "pilot"  code.  Never t  lie  1  ess  ,  the  unique 
. i r nr  i  1  i ;  i s  of  ESP  have  b"en  shown  to  provide  substantial  advantages  over 
'“us  approaches  to  the  store-flutter  problem.  Therefore,  this  user's 
i  .  ;1*.  was  prepared  to  permit  early  utilization  of  ESP  on  practical  problems. 

,  \ 
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2  -  INTRODUCTION 


During  the  development  of  the  initial  version  of  the  Flutter  And  STrength 
Optimization  Program  (FASTOP),  Reference  1,  under  contract  F3361 5-72-C-l 101 
from  the  Air  Force  Flight  Dynamics  Laboratory,  Mr.  Keith  Wilkinson,  the 
project  engineer  on  that  contract,  recognized  that  much  of  the  technology 
being  used  for  minimum-weight  structural  resizing  in  FASTOP  also  had  the 
potential  for  substantially  reducing  the  time  and  cost  required  to  determine 
vhi.h  combinations  ot  wing-mounted  external  stores  would  result  in  the  lowest 
■  :r.  r,P  r  flutter  speeds. 

Subsequent !v,  under  contract  N00019-76-C-0160  from  the  Naval  Air  Systems 
'■'mir.irul,  as  well  as  a  complementary  Grumman  Independent  Research  and  Develop- 
project,  a  search  algorithm  for  wing/store  flutter  was  developed,  refined 
■so  tested  bv  modifying  and  expanding  the  FASTOP  code  (see  References  2  and 
V.’Ven  this  development  effort  led  to  a  pilot  program  that  exhibited  both 
g->cd  *•,  liability  (absence  of  search  failure)  and  good  convergence  character- 
'•  ••  :  it  ■>  work  was  continued  under  a  second  NASC  contract,  N000 1 9-79-0-0062 ,  to 
>>d-l  b’jt'ju's  desirable  for  practical  applications  and  to  demonstrate  the  new 
i.x  t  ,■  r  !;,a  1  -  Stores  Program  (FSP)  on  a  representative  attack  aircraft  and  its 
■  u  i  ceil  store  inventory  (set  Reference  A).  The  project  engineer  on  these 
•:ud  t  s  was  Mr.  Richard  Ghipman. 

With  the  performance  and  the  advantages  of  the  store-search  procedure 
i v • ; , g  hfc.t  confirmed ,  utilization  of  the  procedure  on  current  aircraft 

'eels  hr- are  desirable.  To  permit  early  availability  of  FSP  to  practicing 
■  analysts,  ami  to  increase  the  applicability  of  the  program  to  modern 
. .  i  *  craft  with  thinner,  more  flexible  wings,  the  following  additional 
.  -  taken  under  contracts  NOOO  I  9-8  1  -C-0395  and  NOIhi  1 9-84-G-01  2  1 : 

’  Wvoni  technical  enhancements  were  introduced,  including: 

An  Increase  from  20  to  40  in  the  maximum  number  of  modes  that  can 
be  calculated  in  a  vibration  analysis  and  used  in  a  flutter 
anal v s i s 

i':  -iv  is  ion  lor  including  rigid-body  modes 


o  Logic  for  eliminating  modes  from  a  flutter  analysis  based  on 
ratios  of  generalized  forces  to  generalized  masses 
o  An  increase  from  6  to  15  in  the  maximum  number  of  reduced  veloc¬ 
ities  at  which  generalized  aerodynamic  forces  are  calculated  and 
later  used  as  a  base  for  interpolation. 

')  Logic  was  added  to  by-pass  data  required  for  store-search  runs  if 
onlv  a  conventional  flutter  analysis  is  desired. 

0  The  pilot  code,  previously  available  only  as  an  IBM  version,  was 
converted  to  a  CDC  version  compatible  with  the  Central  Computer 
System  at  the  Naval  Air  Development  Center.  This  included  a  sub¬ 
stantial  reduction  in  computer  central-memory  requirements. 

41  A  capability  was  introduced  for  reading  the  dynamics-model  input 
matrices  either  directly  from  NASTRAN  output  files  or  from  card 
Images . 

’■  This  user's  manual  was  prepared. 

•.void  unnecessary  rfupl  ica  t  ion ,  the  content  of  this  report  complements 
ir. tended  to  be  used  in  conjunction  with  prior  FASTOP  and  ESP  reports, 
ter  ip t ions  of  input  data  common  to  FASTOP  and  ESP,  the  reader  is  re¬ 
in  Volume  II  of  either  Reference  1  or  Reference  5.  Theoretical  docu- 
i  v.  of  die  vibration-  and  flutter-analysis  methods  is  contained  in 
i  af  Reference  1 .  Theoretical  documentation  of  the  store-search 
:  is  contained  in  Reference  4. 

version  of  ESP  described  in  this  report  was  developed  for  execution 
-miputing  equipment  having  extended-memory  capability,  and  most  of  the 
.  oti  contained  herein  should  be  applicable  to  any  such  installation. 

■  ,  the  ■’ufo;  mat  ion  on  control  cards  in  Section  4  is  written  specifically 
;  equipment  and  NOS  operation  system  at  the  NATO  Central  Computer 


3  -  PROGRAM  OVERVIEW 


The  External-Stores  Program  (ESP)  performs  a  gradient-directed  numerical 
search  to  determine  those  combinations  of  external  stores  that  are  most 
flutter-critical  for  a  given  aircraft  design  and  flight  condition.  The  search 
is  performed  in  a  parameter  space  defined  by  a  user-specified  range  of  possi¬ 
ble  store  properties  (including  both  mass  properties  and  on-diagonal  pylon 
r  f  h  i  1  :  t  i  es  ,  if  desired)  at  each  store  station.  A  search  step  consists  of  a 
A-  •  i*  ion-mode  analysis,  a  flutter  analysis,  a  calculation  of  the  derivatives 
of  the  t  hitter  speed  with  respect  to  the  various  store  parameters,  the 
deti  i mi  tin r ion  of  a  new  set  of  store  parameters,  and,  to  start  the  next  search 
step,  a  nodi  f  leaf,  ion  of  the  system  mass  and  flexibility  matrices  based  on 

•  hrse  new  parameters.  The  search  is  started  from  a  user-specified  set  of 

.--ue  parameters ,  and  is  continued  until  either  a  local  minimum  flutter  speed 
-  ;;  led  user-specified  number  of  search  steps  has  been  performed.  If 

am  er  uce  to  a  minimum  has  not  been  achieved  after  the  specified  number  of 
■ -u  s  ,  a  file  containing  recent  values  of  key  search  parameters  is  saved  to 

■  s'  in  continuing  the  search  in  a  new  job  submission  if  desired.  To 

•  • ;  ,*  us-r  assessment  of  the  progress  made  during  a  search,  an  auxiliary 

prepared,  separate  from  the  primary  output  to  he  printed,  in 
! :  ;  s  summary  is  given  of  kov  parameters  at  each  search  step. 

’  !  i’.RA  f  ION-ANALYSIS  MODULE 

•  v'. lr.it  ton-analysis  procedures  in  ESP  are  basically  the  same  as  those 

latter  Ant!  STrength  Optimization  Program  (FASTOP)  from  which  it  was 
u  !  .  however,  the  means  of  providing  most  of  the  dynamics-model  input 

’  ■ ’  has  he  on  changed  to  facilitate  the  transfer  of  this  information 
■-cam  programs  other  than  the  Strength  Optimization  Program  (SOP)  in 
'-'p-  ri  f :  rally,  there  is  now  a  direct  interface  to  ESP  from  both  the 
1  i  ■  Lhe  MacNea 1 -Schwend 1 ei  Corporation  (MSC)  versions  of  NASTRAN  for 
i  1  ii  !  t”  matrix,  two  matrices  containing  the  required  mass  data,  and  a 
•  d i sp Incements  in  the  dynamic  degrees  of  freedom  due  to  unit  rigid- 
o  ,  -  1  a  ■  orient  s .  A  description  of  the  steps  required  in  a  COSMTC  or  MSC 
'  '0  1  u>  *<}  obtain  these  matrices  in  an  ESP-compatible  form  is  given  in 

■  :  -  A  ''tie  dvT'«un  i  cs-rr-ulel  input  matrices  mav  also  lie  read  into  F.SP  as 


:.'i  -image  files,  which  could  be  created  either  manually  or  by  a  user-provided 
-conversion  program  that  is  run  following  the  execution  of  another 
•ream  program. 

other  ESP  vibration-analysis-module  changes  from  FASToP  include  an 
. ease  in  the  maximum  number  of  dynamics-model  degrees  of  freedom  (220  vs. 

>)  and  vibration  modes  (40  vs.  20).  Also,  as  a  user  option,  rigid-body  mode 
tt'-'s  can  be  appended  to  the  structural  vibration  modes  for  inclusion  in  the 
■  •(■  rent  flutter  analysis;  the  maximum  total  number  of  modes  is  also  40.  To 
,.n:r  limited  degree  of  user  control  over  rigid-body  dynamic  characteristics 
if  cannot  be  modelLed  well  or  at  all  with  the  present  capabilities  of  ESP 
• ,  aerodynamics  of  nonconvent  ional  fuselage  geometries,  or  flight-control  - 
-  -  ci. vets),  norite  ro  values  may  be  specified  for  the  zero-airspeed  fre- 
■  or  the  rigid-body  modes.  This  capability  will  provide  maximum 
’  i  1  wh u  ;ca  in  conjunction  with  the  capability  in  the  flutter-analysis 
■  •'  ;■  r  ov  ious  Ly  available  in  FASTOP)  for  specifying  zero-airspeed  modal 

o.i.  '  liiies.  It  Is  suggested  that  the  input  values  be  selected  so  that, 

■■■■ ined  with  the  effects  of  the  program-calculated  aerodynamics,  the 
:  : ■! -body  characteristics  at  velocities  close  to  the  flutter  speed  are 
A  discussion  ot  the  implementation  of  the  rigid-body  modal 
:  s  t  y  in  the  ESP  vibration-analysis  module  is  contained  in  Appendix  B. 

1  'EK-ANA(,iStS  MODUI.E 

; ter -analysis  procedures  in  ESP  also  are  basically  the  same  as 
■..\vrnp,  nut  again  some  significant  changes  have  been  introduced. 

’  :  t  lie  Inc  t  ease  in  the  maximum  number  of  modes  in  the  vihration- 
the  corresponding  number  for  the  flutter-analysis  module  also 
'  '  c  •  -  .•  i  torn  2  If  to  40.  In  addition,  to  address  the  much  broader 
••'luce.'.:  velocities  that  must  he  used  when  rigid-body  modes  are 
’  "  f  iut  c  -:i  analysis,  a  change  from  6  to  15  was  made  in  the  maximum 

••••un.-.-ul  velocities  tor  which  generalized  aerodynamic  forces  may  be 
•  :  • i ect ly  and  biter  used  as  a  base  for  interpolation.  Related  to 

’  1  th--  general  i  zed-force  interpolation  accuracy  test  described  in 

Volume  n,  page-:  88  through  9’,  has  been  deleted,  and  the  actual 


number  of  reference  reduced  velocities  to  be  used  in  a  particular  analysis  is 
now  a  user-specified  quantity. 

Another  flutter-analysis-module  modification  related  to  the  capability  to 
Include  rigid-body  modes  was  a  change  to  a  two-zone  general i zed-foree-i nterpo- 
lut  ion  method  (one  zone  for  high  reduced  velocities  and  one  for  low  reduced 
velocities)  in  the  flutter-solution  procedures.  With  this  approach,  no  values 
•t  i :,e  independent  or  dependent  variable  approach  infinity  for  either  very 
del;  ■  very  low  reduced  velocities.  This  change  is  transparent  to  the  user. 

A  n.-w  feature  introduced  into  the  flutter-analysis  module  is  a  capability 
.  mt  op.  it  ly  e  1  iir.i  nut  i  ng  from  a  p-k  flutter  solution  any  modes  that 
n...  !  d  nor  ■  :  on  i  f  leant  ly  affect  the  predicted  flutter  characteristics.  A  mode 
e  i  'P'i  .Kited  vhen  the  ratio  of  the  modulus  of  the  on-diagonal  genera  H  zed- 
. e  nn  in  the  value  of  the  corresponding  generalized-mass  term,  is  below  a 
■  -spec  i  f  ■’  e .:!  value  of  that  ratio  for  each  of  three  airplane  forward  ve lo¬ 
ci  i  i  e s .  ''luce  the  computing  time  used  in  the  p-k  solution  is  approximately 
n  i  it  i'in.il  to  the  cube  of  the  number  of  modes,  and,  for  problems  with  many 
,  he  p-k  Solution  time  is  a  large  percentage  of  the  total  time,  utiliz- 
n,;  :!  i  feature  can  significantly  reduce  computing  expenditures. 

Arc’ll  :h)i:i  f.E 

■  h  ■  search  objectives  of  ESP  and  FASTOF  are  verv  different,  much  of 
-.lie  F  SP  flutter-optimization  module  is  new.  Flutter-velocity 
t  .  -  ;  ji'..  .;alculated  tor  store  mass,  inertia,  and  center-of-gravi  ty 

•  '  ’••is  ,r\i  pylon  flexibilities  for  up  to  five  store  stations,  and  a 

•  *  •  .''too  niun»T  leal  search  technique  is  used  to  determine  the  loca- 
- 1  ore  -pa  rasne  t  «.•  r  /  py  1  on-f  1  ex  ib  i  lity  space  at  which  the  flutter 
,  ’  •  .•  lr  •  m  :i- .  A  detailed  description  of  the  ESP  search  procedure  and 

•  .lire '-at  ions  is  given  in  Reference  4,  Sections  3  through  5. 

•  '  :  ■  ’  ‘  I  V  OPTION 

1  ■  '  •  •  u  ’  .  b-ei  tivi  o*  developing  KSP  was  to  provide  an  automated 

•  v  dui i  ' '  i  determining  flutter-critical  combinations  of  external 


Ills  program,  like  FASTOP ,  also  may  be  used  in  a  conventional  analysis 
t sired.  By  setring  a  single  clue,  program  steps  involving  a  search 
seed ,  and  input  data  required  for  searches  can  he  omitted, 

;  RAM  St  ATI'S 

tri^e'ii  at  FSI’  documented  in  this  report,  although  erdiarced  in 
e i"  ,  c  s  trout  the  version  that  existed  at  the  conclusion  of  the  work 
:  n  reference  A,  must  still  he  considered  a  "pilot"  code.  The  program 
res  not  been  optimized  for  minimum  computing  time,  and,  for  the  most 
1 h  use  options  actually  used  for  the  demonstration  described  in 
.  <.■■_>  seen  checked .  For  the  options  that  have  not  been  checked, 

:  that  effect  have  been  included  as  footnotes  to  the  corresponding 

:  -.c  ,  i  u  1 1  <  uis  in  Subsection  5.2, 


it  is  noted  chat  the  data  descript  ion  for  FASTOP-  1  given  in  Reference  S , 
.,e  :i,  pages  224-307,  is  identical  to  the  description  in  Reference  1 
,•;>(  for  KI.UE(8)  in  Item  6  of  the  main-program  data,  which,  as  specified 
,  i i  dummy  variable  in  FSP.  Therefore,  the  following  text  may  be  used 
,  ■  i;  ;  line  t  i  or.  with  either  of  these  references.  To  assist  in  cross¬ 
es  c  i  r.g  ,  the  page  numbers  in  References  1  and  5  corresponding  to  each 
i  -.roup  of  data  have  been  listed  in  the  following  subsections  immediately 
List  :  r  titles. 


1  •  .*  L  a  Entered  Via  Main  Program 

'See  also  Reference  1,  Volume  1 ! ,  pages  195-201,  or  Reference  5, 
,'o  1  jmt.  IT,  pages  2  2  ’  -  2  2 <4 .  I 

5  Identical  to  FASTOP . 


k:  PP,  :) 

1  dent  teal  to 

FASTOP. 

is  l.I7  K  ( ■?< ) 

•  h  ■  k  ( ; '  «  o 

i a >  not  enter 
perform  nr.aly 

store-search  module,  i.e., 
sis  only. 

7 

Periorm  .-tore 
i  and  PlPFfa! 

search.  Requires  that  KM'K(3) 

:  4 . 

Ml  i:  i  8 )  =  0 

’  '..n."".  variab! 

e  in  FSP. 

VI  : '  (•  ■'  a  ) 

identical  to 

FASTOP. 

K  i (  .;  }  : 

pi  "F  <  :>•)  -■  " 

f’lnwrv  '’or  ,0,1 

e  in  FSP. 

[ 1 .  _  J  i 

\  »  /'  lh.:  ; 

n  FSP. 
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a  i  i'.1  if.  ISP. 
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in  P-P. 
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i  The  first  data  group  on  each  card  should  be  nonblank;  one  or  both  of 
the  other  two  tnav  he  blank. 

v  few  additional  comments  on  the  dynamics-model  matrices  follow.  First, 
noted  that  the  modal-interpolation  procedure  in  the  flutter-analysis 
•  requires  that  the  input  modal  data  he  specified  at  points  along  a  set 
,.nwi  se—'r  United  lines  (see  Reference  1,  Volume  I,  page  P')  .  Therefore , 
r  ‘tiics  idealization  used  to  create  the  input  matrices  to  the 
: on-sral ys is  module  should  be  defined  with  this  requirement  in  mind. 

•frC'-.nd  the  vibration-analysis  module  in  ESP  cannot  accommodate  a 
ar  mass  matrix.  Therefore,  the  masses  at  all  dynamics-model  grid  points 
>i  'zero,  and,  at  these  grid  points  for  which,  rotation  degrees  of 
u  are  used,  the  corresponding  moments  of  inertia  also  must  he  nonzero. 

when  NASTRAN  if  used  tut  define  the  FSP  dynamics-model  matrices, 
•rrU  es  he  passed  (using  the  procedures  described  in  Appendix  A)  are 
•u.r : esronding  to  the  ’analysis"  (or  "solution”)  set.  Therefore,  the 
'  :.  -u-  ;i  1  tzatiem  used  must  he  such  that  the  store  degrees  of  freedom  are 

r.denr  nr. rd inatcr. . 

I'.IITT-DATA  FII.F 

:  i  .  •  he  input  dr-tn  to  fSP  is  to  a  considerable  extent  similar  to  that 

'■■i!  for  U.STU’’  in  Reference  1,  Volume  II,  pages  195-279,  the  primary 
is  in  the  text  below  is  on  data  which  is  new  for  ESP.  To  avoid  unneces- 
'■  ,  1  :  r  .t  e  blocks  of  descriptive  information  for  FASTOP  which  are 

•  :  t  r  ,  or  ate  not  used  in  ESP ,  are  simply  described  as  such. 

•i  !  '  i  at .  rn  ss-re  fo  renc  ;  rig  ,  the  data-item  numbers  in  Reference  1 
■  •» :  t-.-ui  below.  Wlier-  a  data  item  that  was  previously  used  in  FASTOf 
•o  ' if  i.  "I  for  use  in  K  Si' ,  the  associated  variable  i  ist  and  description 
Id  c  nsidertd  to  replace  the  information  in  Reference  1. 

.  - 1 1  •  •  •  i  hat  have  been  retained  unchanged  in  FSP,  a  new 

■  1  ,  i  t  ;  [■  has  been  given  t  ■  provide  further  cl  aril  i  cat  ion.  For  da  t  a 
•  v  t  ot  1  .  i'.'w  i  1.  em  numhe’-  ;  have  hern  created  Iv  appending  'otters 

:  b  i :  I  A  fO. 1 1  !  •  eir  numbers. 


ow  the  dynamic  mass  matrix  in  input  file  TAPE25.  In  the  single  input  file 
for  COSMIC-NASTRAN  matrices  (TAPE20) ,  the  following  order  must  be  adhered 
flexibility  matrix,  mass  matrix,  rigid-body-displacement  matrix,  plug 
matrix.  Additional  specific  instructions  for  writing  the  NASTRAN  files 
contained  in  Appendix  A. 

When  the  dynamics-model  matrices  are  provided  in  card-image  form,  the 
owing  data-preparat ion  instructions  may  he  followed  for  all  three 

i  e  e  s  . 

.  The  Fortran  format  should  be  (3(214  ,F.l  S . 5, IX) )  .  Each  group  ot  three 
numbers  in  this  format  consists  of  a  row  number,  a  column  number,  and 
■  in  element  value. 

,  A'1  matrix  elements  should  be  provided. 

rhe  elements  should  he  in  consecutive  older  by  row. 

The  first  element  of  each  row  should  constitute  the  first  data  group 
on  a  new  card. 

Ka;  h  card,  except  that  containing  the  last  element  in  each  row,  should 
have  a  full  complement  of  three  data  groups. 

The  cards  containing  the  matrix  elements  should  be  followed  by  a  blank 

c  it  rd . 

i  resiled,  die  instructions  just  listed  may  be  relaxed  for  individual 

■ as  follows: 

•  ■. 1  b  •'  ]  i  r  v  matrix: 

Tt..  row  and  column  numbers  are  optional;  i.e.,  the  format  may  be 
!  V8X.F1  S.  S,  IX)  1  . 

final  blank  card  is  not  required. 

•  "•  r  V  '■  X  l 

i  '■  mnii'M  matrix  elements  must  be  provided. 

■  I  cl’..  “  r  i  r « ■  t  ■  data  group::  on  a  card  may  be  blank,  but  each  card 

*>>••■  :t  tl.e  last  (blank)  card  should  have  at  least  one  data  group, 
i  1  rr-t  element  specified  tor  each  row  should  he  on  a  new  card. 

■i  -  h  >■  I  v-d  i  so  1  acemei.  t  matrix: 

n,  :  re:, i  matrix  elements  must  he  provided . 

*  ■;  r.o  >  !  ement-order  ing  renti  i  remunt . 


till-  inboard  direction,  and  positive  is  associated  with  upward  vertical 
i  placement  of  the  left  side.  These  modifications  produce  a  sign  change  in 
I  I  off-diagonal  terms  in  the  sixth  row  and  column  of  Eqs.  (4-2)  through  (4-5) 
f  Reference  4.  The  degrees  of  freedom  that  should  adhere  to  the  coordinate 
vstem  in  Figure  5-1  are  primarily  those  that  are  used  at  the  pvlon/store 
:  uchment  points.  However,  for  a  wing  represented  structurally  as  a  beam, 
be  decrees  of  freedom  that  are  used  for  modal  interpolation  in  the  flutter- 
.lysis  module  should  be  defined  such  that  nose-up  rotation  in  pitch  is 

■  isfdert i  positive  when  downward  vertical  displacement  is  positive. 

■hen  a  store-parameter  search  is  being  performed,  certain  elements  in  the 

■  .  I  j 1 u  v  and  mass  matrices  that  are  supplied  to  ESP  may  be  arbitrary.  For 
!,.>  tie..  Ibilify  matrix,  the  on-diagonal  elements  associated  with  the  store 

cs  of  freedom  are  arbitrary,  and,  for  the  mass  matrix,  all  store-degree- 
•  •  r  . . ion.  elements  are  arbitrary.  Before  the  first  vibration  analysis  is 


Figure  5-1  -  Store  nvnamir  Coordinate  System 


5  -  INPUT-DATA  DESCRIPTION 


5 . !  -  DYNAMICS-MODEL  DATA  FILES 

As  indicated  in  the  previous  section,  there  are  four  dynamics-mode] 
matrices  that  must  he  provided  to  ESP  from  an  external  source: 

•  ->  A  flexibility  matrix.  This  has  been  called  KLLI  in  the  illustrative 
NASTRAN  DMA?  alter  statements  shown  in  Appendix  A  in  Figures  A-2  and 
A- 4. 

o  A  mass  matrix  associated  with  the  dynamic  degrees  of  freedom  other 

than  those  that  are  assumed  to  be  fixed  when  computing  the  flexibility 
matrix.  For  brevity,  this  matrix  will  hereinafter  be  referred  to 
slir.piy  as  the  dynamic  mass  matrix  or  mass  matrix.  The  NASTRAN  name 
used  for  this  matrix  in  the  DMAP  alter  statements  in  Figures  A-2  and 
A-h  is  MLLW. 

e  A  separate  mass  matrix  for  the  degrees  of  freedom  at  the  assumed 

frec-body  support  points.  Following  the  terminology  in  Reference  1, 
Volume  f,  Section  ?,  pages  48  and  49,  this  will  hereinafter  be 
referred  to  as  the  "plug"  mass  matrix.  The  name  MRRW  has  been  used 
for  this  matrix  in  the  illustrative  DMAP  alter  statements  cited  above. 
4  matrix  containing  displacements  in  the  dynamic  degrees  of  freedom 
die  to  unit  rigid-body  displacements  at  the  free-hody  support  point  or 
plug.  This  matrix  is  called  DM  in  both  the  COSMIC  and  MSC  versions  of 
NASTRAN  Rigid  Format  5. 

'■■■  '  r  matrix,  as  read  In  ESP,  is  arranged  such  that  the  displacements  in 

it;  various  dynamic  degrees  of  freedom  due  to  each  unit  rigid-body 
!:  a  i  -  cense:.:  occupy  one  row  of  the  matrix.  The  rows  should  he  ordered  such 
:  hr:  i  'iji1  oi  ly  translation  modes  precede  rigid-body  rotation  modes.  The 

i i  'his  matrix  is  the  >  submatrix  in  Fq .  (7.15)  of  the  reference 

'•u  cam!  r  s -model  matrices  should  apply  to  a  half-airplane,  and,  for 
: ’  :n  r"V  degrees  of  freedom,  should  be  developed  following  the  coordinate 
.  1  i  r  vn  in  F-  gure  S-]  .  This  coordinate  system  is  a  variation  of  that 
r  '  igure  4-1  of  Ref erer.ce  4:  To  permit  utilization  of  matrices  from 
'A:  ,  wo:  ;’i  ir<  based  in  a  right-hand  coordinate  system,  Y  ir.  now  positive 


5-1 


As  a  final  comment  on  the  illustrative  control-statement  sequences.  It 
I'te.i  chat  the  implied  judgment  concerning  which  files  are  to  be  considered 
i  i rec t -access  and  which  are  indirect-access  should  be  used  lor  initial 
nii dance  only.  Also,  where  there  is  a  significant  possibility  that  a  file 
ould  he  in  either  category,  depending  upon  its  contents  for  a  particular 
implication,  the  direct-access  option  has  been  selected.  This  is  especially 
ipplicuble  to  TAPF24,  the  mass-matrix  file.  if  this  matrix  is  obtained  from 
'■  NASTRAN ,  it  can  probably  be  considered  as  indirect-access,  whereas,  if  i 
s  in  .aid-image  form,  it  might  have  to  be  direct-access. 


TAPF08  is  a  composite  of  up  to  15  files  which  contain  aerodynamic- 
influence-coefficient  matrices  (see  note  1  on  page  4-5)  calculated  at  a 
specific  Mach  number  for  the  reference  reduced  velocities.  In  the 
illustrative  control-statement  sequence  for  an  analysis  run,  shown  in  Figure 
4-1,  Arc's  are  being  created  and  saved,  and,  in  the  two  search  runs,  shown  in 
Figures  4-2  and  4-3,  these  AIC's  are  being  used.  The  creation  of  AIC's  in  a 
separate  initial  analysis  run,  rather  than  as  part  of  a  search  run,  is 
general !v  recommended,  since,  in  a  search  run,  the  AIC's  would  be  unneces- 
,.ii  !  lv  recreated  during  each  step  in  the  search.  As  a  further  measure  to 
<  .••user,  e  computing  time,  it  is  suggested  that  the  initial  analysis  run  utilize 

•  1,,-  V  method  for  the  flutter  solution  rather  than  the  p-k  method.  A  quantita- 
>  •  o  : uc.  1  cat ! on  of  the  savings  to  be  realized  from  these  guidelines  can  be 
■ot-iined  from  Table  4-1. 

i he  rifth  rossible  input  file,  read  via  TAPE47 ,  is  used  only  when  restart 
fug  (see  Figure  4-3).  This  file  would  have  been  created  via  TAPE48  as  an 
output  file  fr  >m  a  prior  search  or  restart  run.  Additional  information  on 
•ui-.i  iile  is  contained  in  Subsection  6.3,  beginning  on  page  6-24. 

'  j.  to  four  output  files  may  he  created  and  saved  via  ESP.  Tn  addition  to 
f  he  A I  •  •  file  and  the  restart  file  which  have  already  been  discussed,  there  is 
!  ■  ,  fill'  for  obtaining  plots  of  the  flutter-analysis  results,  and  a  file 

i  i  •<  i  ng  a  summary  of  key  parameters  at  each  step  of  a  search.  The  latter 
1  ;  Trained  via  1APF40  an  Illustrated  in  Figures  4-2  and  4-3.  In  most  cases, 

.  qei'-V  interactive  review  of  this  file  will  permit  a  search  to  be  restarted 
v ' r h  confidence  If  required. 

lie  plot  iile,  which  may  he  obtained  from  any  search  or  analysis  run,  is 

*  • to-  tin  TaPK 60 .  in  a  search  run,  the  plot  file  contains  the  plotting 

a  i  p  i  it  a!'.  I  hitter  analyses  performed  during  the  search.  The  conversion 

(.ho  plot  file  to  actual  hard-''opv  plots  using  an  NADC  CialComp  plotter  is 
or r  'me : f shed  via  a  manually  completed  request  form.  Therefore,  this  step 
t  c<  u  :  .-e'.  *  he  assistance  of  an  on-site  NADC  person. 


I  tie  two  parameters  on  the  ACCOUNT  cards  are  a  user  number  and  a  password. 
Lei'!  should  be  requested  as  described  on  page  1-22  oi  Reference  6. 

An  absolute  file  of  the  ESP  program  is  accessed  via  the  first  ATTACH 
td.  The  four  digits  in  the  file  name  designate  the  version  of  the  program 
terms  of  the  year  and  the  month  in  which  it  was  created.  The  six-digit 
.  omi  parameter  on  the  program  ATTACH  card  designates  the  user  number  under 
ich  the  program  absolute  file  is  stored.  The  third  parameter  is  the 
i-  'ciated  password.  Navy  personnel  wanting  to  run  ESP  at  NADC  may  obtain  the 
st  two  parameters  by  contacting  Mr.  Robert  Richey,  Code  Ci04  1  2 ,  Naval  Air 
,e  opment  Center,  Warminster,  PA  18974,  (21b)  441-1944. 


Input  data  to  ESP  consists  of  a  user-prepared  card-image  file  read  via 
'E!6,  plus  up  to  five  additional  files  depending  on  the  options  specified  in 
i'PE'b  file  From  one  to  three  of  these  files  -  read  via  TAPE20,  TAPE25, 
TAPEl’b  -  contain  four  matrices  which  together  constitute  most  of  the  data 
..id  to  define  the  analytical  model  of  the  airplane  dynamics.  These  are: 
U-xibi t i ty  matrix,  a  mass  matrix  associated  with  the  dynamic  degrees  of 
vdimi  other  than  those  that  are  assumed  to  be  fixed  when  computing  the 
j:.!1!'.  i;ty  matrix;  a  separate  mass  matrix  for  the  degrees  of  freedom  at  the 
ur.r.cd  tree-body  support  point;  and  a  matrix  containing  displacements  in  the 
'•aii'  degrees  of  freedom  due  to  unit  rigid-body  displacements  at  the  free- 
I,  support  point.  Each  of  these  matrices  may  be  read  into  ESP  In  a  form 
i.nt  with  output  from  NAS  TRAN  (either  the  COSMIC  or  the  MacNeal- 
,d!  er  Corporation  (MSC)  version),  or  they  may  be  read  as  card-image 
Matrices  obtained  from  COSMIC  NASTRAN  are  all  read  via  TAPF.20. 

• -CAS IRAN  matrices  or  card-image  matrices  are  read  as  follows:  the  flexi- 
r  r  I  y.  via  TAIE20,  the  two  mass  matrices  via  TAPE2S ,  and  the  rigid- 
••  displacement  matrix  via  TAPE26.  In  the  illustrative  control-statement 
...  i  .  s  shown  in  Figure  4- i  through  4-3,  it  has  been  assumed  that  all 
•  1  re  •’  are  coming  from  COSMIC  NASTRAN,  and  therefore  only  TAPE20  is  active. 

b- '  •  s'a -TRAN  or  card-image  matrices,  the  statements  that  have  been  corn¬ 
s'.1  "it  with  asterisks  would  be  activated  in  place  of  the  currently  active 
!  •. ra rerneiu  . 


Table  4-1 


Typical  Central-Processor  Times 
tor  Major  Components  of  an  ESP  Execution  Run 


Kune t ion 


:  brat  i  on  .it.n  1  ys  is 

Compute  and  save  double-lattice 
aerodvn.'imic  Influence  coef- 
!  feients  (AFfi's)  and  compute 
generalized  forces 

; 

previously  saved  AIC's'  to 
'•'iRp.it e  generalized  forces 

V -me *•  hod  flutter  solution 


j.  -k  met!  rd  flutter  solution 
p-K  method  flutter  solution 

nivtv  j  ,1  tit  ter  solution 


•  •  r !  vat  i  .»*>  calculation  and 
■  ■  ' : ch  slip 


Cyber  760 

Problem  Size  CP  Time 

(dec.  secs. 


200 

degrees 

of  freedom 

90 

40 

flexible  modes 

60 

panels 

60 

15 

reduced 

velocities 

8 

modes 

60 

pane  1 s 

6 

15 

reduced 

velocities 

8 

modes 

20 

reduced 

velocities 

1 

8 

modes 

40 

velocities^ 

70 

8 

modes 

40 

.  .  4  3 

velocities 

6000 

40 

modes 

12 

veloc i ti es^ 

3400 

40 

modes 

200 

degrees 

of  freedom 

1 

8 

modes 

6 

store  parameters 

•  niri  saved  are  not  AIC's  but  intermediate  results  of 
i " r inn  procedure  i see  Reference  7,  pages  34-35). 

■'"'iia’ ions  of  derivatives  of  generalized  forces  with 

i  e  li.er-il  f  reeu.  my  . 

:  .  r  st  ei  it  search . 


ESPRUNR,TXXXX,CM300000»EC70»STPOO. 

ACCOUNT  , UUUUUU , PPPPPPP . 

RFL , EC  =  70 • 

ATTACHf ABS=ESP8501/UN=NNNNNN»PW=WWWUWWU. 
ArTACH,TAPE08=AAAAAAA. 

GET » TAPE  16  =  DDBDliDD. 

ATTACH,TAPE20=CCCCCCC. 

*ATTACH,TAPE20=FFFFFFF. 

*AT  rACH,TAPE25  =  MMMMMMM. 

*GET  ,TAPE26=TTTTTTT  » 

GET  ,  TAPE47=IIIIIII < 

OEFINE»T APE60=G6GG6GG » 

COPYBF  , TAPE08,FL0801A. 

CUP YBF »  TAPE08 » FL0802A . 

COPYBF , TAPE08,FL0803A. 

COP YBF , TAPE08 ,  FL0804A ♦ 

COPYBF  , TAPE08  >  FL0805A  « 

COPYBF  ,  T  APE08 » FL0806A . 

COPYBF, TAPE08,FL0807A. 

COP YBF ,  T  APE08 , FL0808A ♦ 

COPYBF »TAPE08,FL0809A* 

COPYBF, TAPE0B,FL0810A. 

COF  YBF , TAPE08,FL0811A. 

COPYBF , TAPE08 » FL0812A ♦ 

COPYBF ,TAPE08,FL0813A. 

COF YBF  >TAPE08»FL0814A, 

COPYBF, TAPE08,FL0815A. 

REIJIND«FL0801A»FL0802A,FL0803A»FL0804A,FL0805A,FL0806A 

REUiMD,FL0807A,FL0308A,FL080?A,FL0810A,FL0811A,FL0812A 

REVINO»FL0813A»FL0814A»FL0815A* 

ABS » 

SAVE , TAPE40-SSSSSSS . 

SA'.'E,  rAPE48=0000000  . 


e  -  Typical  Control-Card  Sequence  for  Search-Restart  Run  Using 

Previously  Saved  Aerodynamic  Influence  Coefficients. 


ESPRUNS»TXXXX»CM300000»EC70fSTPOO* 

ACCOUNT >UUUUUU»PPPPPPP. 

RFL  f  EC  =  70 • 

ATTACHrABS=ESP8501/UN=NNNNNN.PW=WUWUWWW. 

ATTACH »TAPE08=AAAAAAA. 

GET  » T  APE  1 6  =  DDDDDDD . 

ATTACH»TAPE20=CCCCCCC. 

* ATT ACH r TAPE20  =  FFFFFFF  < 

*ATTACH>T  APE25=MMMMMMM . 

♦  GET » TAPE26  =  TTTTTTT . 

DEFINE»TAPE60fGGGGGGG. 

COPYBF rTAPE08»FL0801A. 

COPYBF  »TAPE08>FL0802A. 

COPYBF* TAPE08*FL0803A. 

COPYBF » T APE08 » FL0804A * 

COP YBF  *TAPE08»FL0805A» 

COPYBF  r T APE08 » FL0806A . 

COP YBF  *TAPE08tFL0807A« 

COPYBF >TAPE08»FL0808A. 

COPYBF *TAPE08»FL0809A. 

COFYBF  >TAPE08»FL0810A. 

COPYBF  »TAPE08»FL0811A. 

COPYBF  rTAPE08»FL0812A. 

COFYBF ,TAPE08»FL0813A. 

COPYBF  >  TAPE08 » FL0814A . 

C CP YBF , TAPE08.FL0815A. 

REUINDfFL0801A,FL0802A*FL0803A»FL0804A*FL0805A*FL0806A 
REUIND*FL0807Af FL0808A,FL0809Af FL0810A*FL0811A*FL0812A 
REWIN D*FL0813A»FL0814A»FL0815A* 

ABS. 

SAME*TAPE40=SSSSSSS. 

CAME » TAPE48  =  0000000 . 


•-2  -  Typical  Control-Card  Sequence  for  Initial  Search  Run  Using 
Previously  Saved  Aerodynamic  Influence  Coefficients. 


iSPRUNA,TXXXX,CM30Q000»EG70»STP00. 

ACCOUNT, UUUUUU»PPPPPPP. 
i‘;F_,EC  =  70. 

AT  T  ACH , ABS=ESP8501 /UN=NNNNNN  ,PUI  =  WWWUUWW. 

GET  >TAPE1£=DDDDDDD. 

AT  TACH,TAPE20  =  CCCCCCC. 

*ATTACH,TAP£20=FFFFFFF. 

*ATTACH,TAPE25=MMMMHMM. 

*GET , TAPE26=TTTTTTT. 

.  “  F  •  NF.  ,TAFE08=AAAAAAA. 

GtFINE , TAPE60=GGGGGGG. 

ABS  . 

RE^If'!D?FL0801A,FL0802A,FL0803A,FL0804A,FL0805A,FL0806A. 
CCU'MDjF'  0807A,FL0808A,FL0809A»FL0810A,FL081iA,FL0812A. 
KEUIND»FL0dl3A»FL0814A»FL0815A. 

COPYBF  ,  F  L  0  8  0 1 A  ,  T  A  P  E  0  8  . 

COPYBF ,FL0802A,TAPE08» 

CGPYBF,FL0803A,TAPE08. 

CjPYBF ,FL0804A,TAPE08» 

COPYBF , FL0805A , TAPE08  » 

COPYBF ,FL0806A»TAPE08. 

COPVBF ,FL0307A,TAPE08. 

CQC' Y£<r  ,FL0808A,TAPE08» 

COPYBF , FL0809A , TAPE08  * 

P  Y  8  F  ,  F  L  0  8 1 0  A  t  T  A  P  E  0  8  . 

COPYBF ,FL0811A,TAPE08» 
cnPYBF . FL0812A, TAPE08 , 

COPYBF  ,  FL.08 1 3A  ,  TAPE08  » 

COP JBF , FL0814A , TAPE08 . 

'■  CP  'B F  ■  FL0815A,TAPE08. 


-  Tvpicai  Control-Card  Sequence  for  Analysis-Only  Run  in  Which 
Aerodynamic  Influence  Coefficients  Are  Calculated  and  Saved. 


4  -  CONTROL  CARDS 


Execution  ot  ESP  on  the  Central  Computer  System  at  the  Naval  Air  Develop¬ 
ment  Center  can  be  accomplished  using  NOS  control-card  sequences  similar  to 
those  shown  in  Figures  4-1  through  4-3.  The  first  figure  illustrates  an 
analysis-only  submission  in  which  aerodynamic  influence  coefficients  (AIC's) 
arc  to  he  calculated  and  saved.  The  second  figure  illustrates  a  search  run 
utilizing  previously  saved  AIC's,  and  the  third  figure  shows  a  restart  of  a 
search  tun.  The  last  letter  of  the  job  name  on  the  JOB  control  cards  in  these 
liginn  wav  chosen  as  an  identifier  for  each  of  the  three  types  of  runs.  The 
control  cards  for  a  fourth  type  of  submission  -  an  analysis-only  run  using 
previously  saved  AIC's  -  would  be  the  same  as  those  shown  in  Figure  4-2  except 
M  •  i '  the  ^sr  two  cards  may  be  omitted. 

Ihf-  central-processor  time-limits  on  the  JOB  cards  for  submissions  to  the 
a-K  Central  Computer  System  should  be  specified  in  terms  of  decimal  seconds 
on  the  COO  6600.  Appropriate  adjustments  for  job  executions  on  the  Cyber  175 
or  .i  /6(J  are  performed  within  the  system.  (The  additional  control  cards 
t-pec  fried  on  pages  3-2-8  of  Reference  6  should  no  longer  be  used.)  To  assist 
the  use-  in  making  initial  time  estimates,  some  approximate  times  for  the 
'•v.  5 or  components  of  an  ESP  run  are  listed  in  Table  4-1.  The  times  are  given 
.  ■<  t •- o  r.vber  760  seconds  because  this  provides  a  more  immediate  indica¬ 
te'.’  !  thi  practicality  of  a  particular  type  of  run  with  the  computing 

tr  r..r>  ’vatlable.  A  conversion  factor  of  3.7  between  Cyber  760  time  and 
;  :  v-n.  •  s  suggested  on  page  3-2-8  of  Reference  6. 

■  S‘u.(  r id  parameter  on  the  JOB  cards  specifies  the  maximum  amount  of 
r  ■,  itnurv  needed  for  ESP  in  terms  of  octal  words.  This  is  followed  by 
.  :>;u  ,r  enter  t  tor  extended  memory,  given  as  the  maximum  octal  number  of 
!  '■  '■  ,  wort’:- .  'be  actual  request  for  extended  memory  is  accomplished  via  the 
F  ,fp-,n(| ,  '.!!•,  let.  is  the  third  control  card  shown  in  Figures  4-1  through  4-3. 

’•  tl  -I  parameter  or  the  JOB  cards  specifies  execution  under  the  NOS  operat¬ 

ic,  >;  ,;;.eu.  a’,  an  "economy"  charging  rate  as  described  on  page  1-20  of  Refer- 
' .  liie  "PO"  priority  (deterred  execution)  that  results  from  the  use  of 
■It-  ■  Ti'i'tcr  is  required  for  ESP  by  virtue  of  its  centra  1 -memory  usage  and 
i  fictions  given  in  Figure  3-2-1  on  page  3-2-8  of  Reference  6. 


4-1 


t  cm  (i 
cont:  .  ) 


KI.UKn?)  =  0 


=  37 


K l .UK  (  38)  =  0 


-=  3  8 


Can  t  i  1  ever-w  i  n  p,  vibration  analysis  to  be 
performed  in  AVAM. 

Free-free  vibration  analysis  to  bo  perl  armed 
in  AVAM. 

Do  not  include  rip, id-body  modes  in  AVAM 
output  passed  to  AKAM. 

Include  rigid-body  inodes  in  AVAM  output. 

Kil'Ki  18)  is  ignored  by  the  program  if  KiJ'K (  1?) 

0. 


Format  and  comments  on  alternative  approaches  t <->  entering  above 
data  are  identical  to  FASTOP. 


Pita  Filtered  Via  V  i  bra  t  ion-Ana  1  vs  i  s  Module 

(See  also  Reference  ),  Volume  Id,  pages  202— i1 16,  or  Reference  3, 
volume  i I »  pager  ? 30-244.) 


toms 

t 


identical  to  FASTOP. 


KI.HKV (  I  )  —  0 
S!,’  FV(7)  =  0 

F  i.i'KV  (  3  )  =■■  n 

■*  1 

MFEVts)  0 


Fixed  value  in  EST  and  FAS'IOP. 

Do  nor  plot  vibration  modes. 

I  ,  2 

Piot  vibration  modes  on  OaiComp. 

bo  not  list  flexibility  matrix  as  used  in 
e  I  seriv.i  1  ue  so  1  u  t  i  on  . 

i.ist  flexibility  matrix. 

Do  not  list  transformed  mass  matrix  as  used  in 
eigenvalue  solution. 

1  ist  transformed  mass  matrix. 


Dummy  variable  in  FSP. 


r  !  i '  i  V  ■■  n  -  i ! 


I  smii.v  v.i  r  i  ah  1  e  i  n  F.SP  . 


be  not  list  flexibility  matrix  obtained 
from  FAS', KAN  or  card-image  file. 

list  input  flexibility  nm t  r  i  x  . 


(  i  ■  r  ,i  1  I  v  limited  r  i ■  .■m.ib'sis-i 


ci  1  v  runs  . 


item  3 
(cent . ) 


KLUEV(8)  =  0 


Do  not  list  dynamic  mass  matrix  as  obtained 
from  NASTRAN  or  card-image  file. 


=  8  List  input  dynamic  mass  matrix. 

The  rigid-body-displacement  matrix  and  the  plug  mass  matrix  are 
always  printed  as  they  are  read. 

The  format  for  the  above  data,  nd  comments  on  alternative 
approaches  to  entering  these  data,  are  identical  to  FASTOP. 

l  tern  3 A  NDYDOF  Number  of  dynamic  degrees  o!  freedom 

(excluding  plug  degrees  of  freedom).  Maximum 
value  is  220. 

IDYFI.X  =  0  Dynamic  flexibility  matrix  is  provided  as 

card-image  data  according  to  the  format  and 
instructions  given  in  Subsection  5.1. 

=  i  Dynamic  flexibility  matrix  is  provided  in  the 

nonsparse  binary  format  obtained  from  the 
0HTPUT4  routine  in  |he  MacNeal-Schwendler 
version  of  NASTRAN. 

=  2  Dynamic  flexibility  matrix  is  provided  in  the 

format  obtained  from  the  0UTPUT2  routine  in 
the  COSMIC  version  of  NASTRAN. 

l  if  ASS  =  0  Dynamic  mass  matrix  and  plug  mass  matrix  are 

provided  as  card-image  data  according  to  the 
format  and  instructions  given  in  Subsection 
5.1. 

=  1  Dynamic  mass  matrix  is  provided  in  the  sparse 

binary  format  obtained  from  the  0UTPUT4 
routine  in  the  MacNeal-Schwendler  version  of 
NASTRAN.  Plug  mass  matrix  is  proviueijl  in  the 
nonsparse  binary  format  from  OUTPUT4. 

-  1  Dynamic  mass  matrix  and  plug  mass  matrix  are 

provided  in  the  format  obtained  from  the 
0UTPUT2  routine  in  the  COSMIC  version  of 
NASTRAN. 


■pnendix  A  for  additfonal  descriptive  and  illustrative  material 
.Ft  NASTRAN /ESP  interface. 


I  r  e  m  3  A 
(emit.) 


1TRNSF  =  (J 


Matrix  of  displacements  in  dynamic  degrees 
of  freedom  due  to  unit  rigid-body  displace¬ 
ments  is  provided  as  card-image  data  according 
to  the  format  and  instructions  given  in 
Subsect i on  5 . 1 . 

-  1  Rigid-hodv  displacement  matrix  is  provided  in 

the  nonsparse  binary  format  obtained  from  the 
OUTPUT 4  routine  in  the  MacNeal-Schwendler 
v  e  v  s  i  o  n  of  N  A  ST:'  A  N  . 

-  2  Kigid-bodv  displacement  matrix  is  provided  in 

the  format  obtained  l rom  i he  OUTPUT*  routine 
in  the  COSMIC  version  of  N  A  ST  KAN  . 

format  =  (ATS).  Number  ot  cards  is  1. 

Data  are  entered  by  subroutine  I'.AADY. 

r’ASTOr  logic /data  items  not  used  in  ESP. 

I-'gic  Item  -  No  Data 

If  a  store  search  Is  to  be  performed  (K!  IT  O)  -  7),  continue  with 
nstructions  below.  Otherwise  (Kill. (7)  -  ()i  ,  go  to  Ttem  16. 

Repeat  items  ftp  -  81  icr  each  store  station.  Then  add  a  blank 
card.  Maximum  number  o’,  store  stations  is  8. 

Vhc  store-mass  parameters  entered  as  Items  80  -  8!'  will  be  used  in 
the  program  to  replace  the  appropriate  elements  of  the  mass  matrix 
r.  rovit’ed  as  a  separate  data  file  (see  Subsection  5.1).  Similarly, 
"he  f ■ ex i h i ! it"  value;  entered  as  Ttem  8F  will  replace  elements  of 
M,"  f lexibi 1 it\  matrix  provided  as  a  separate  data  file.  The 
■  1  e  1 1.  s  in  these  matrices  -bat  will  be  replaced  are  determined  by 
r  he  Jtgrees  of  freedom  enter  -d  <:s  Item  811. 

■'hen  restart  *  ng  a  st  a*  cli ,  Items  PC  -  Kb  should  be  the  valuer, 
or  respond in?.  to  t'u-  pint  in  the  s  t  ore-pa  rame  ter  space  from  which 
■.he  next  series  of  -e.-n  eh  steps  is  to  begin .  These  are  obtained 
,  car  tie  output  at  the  of  t  lie  previous  ser  ‘ft  a  of  search  steps 

-.ec  Subsection  6.1). 


i  i  l  add l t i opa  ! 


'  rative  material 


tom  fib 


IDSTR 


Store-station  number.  Must  he  one  of  a 
sequence  from  1  to  the  total  number  of  stores. 
Maximum  value  is  5. 


Format  =  (15).  Number  of  cards  Is  1. 


8C 


Data  are  entered  by  subroutine  READY. 

STRWI  Initial  store  weight  in  store-search  run,  lb. 

Format  =  (E15.5).  Number  of  cards  is  1. 


Data  are  entered  by  subroutine  READY. 

tern  8D  STRTI(K),  K=i,3  Initial  store  moments  of  inertia  about  local 

x,  y,  and  z  axes^through  store  center  of 
gravity,  lb. -in. 

Format  =  (3E15.5).  Number  of  cards  is  !. 


BE 


Data  are  entered  bv  subroutine  READY. 


STRRI(K),  K=l,3  Initial  store  center-of -gravi ty  x,  y,  and  z 

offsets  relative  to  store  dynamic  degree- 
of-freedom  location  (nominal  store/pylon 
attachment  point),  positive  for  store 
center-of-gravity  location  forward,  inboard, 
and  below  pylon  attachment  point,  in.  As 
discussed  further  in  Subsection  5.1,  this 
represents  a  change  from  the  convention  used 
in  Reference  4. 


Fo rma  t 


(3E15.5).  Number  of  cards  is  1. 


Data  are  entered  by  subroutine  READY. 

STRSI(K),  K=l,6  Initial  flexibility-matrix  diagonal  elements 

for  translations  in  x,  y,  and  z  directions  and 
rotations  about  x,  y,  and  z  axes  at  nominal 
stcre/pylon  attachment  point,  in. /lb.  or 
rad . /lb . -in . 

These  elements  are  incremented  within  the  program  when  pvlon- 
f'exibility  parameters  are  search  variables.  If  a  flexibility- 
matrix  element  is  entered  here  as  zero,  the  initial  search  value 
of  chat  element  u i 1 1  be  the  one  contained  in  the  flexibility 
matrix  provided  as  a  separate  data  file  (see  Subsection  5.1). 

Format  =  (bElO.l).  Number  of  cards  is  1. 


I 


I 


Data  are  entered  bv  subroutine  READY. 


I  tent  8(.' 


Dummy  variable  in  FSF . 


ISTDOF(K)  =  0, 

K=  1  ,f> 

Format  =  (615).  Nuiuuer  of  cards  is  I. 

Data  are  entered  by  subroutine  RFADY. 

Item  8H  1DYD0F(K) , K= 1 , 6  Dynamic  degree-of-f reedom  numbers  correspond¬ 

ing  to  translations  in  x,  v,  and  ?.  directions 
and  rotations  about  :< ,  v,  and  z  axes  at 
nominal  store/pylnr  attachment  point.  If  an 
airplane  dynamics  model  does  not  include 
degrees  of  freedom  tor  ore  or  more  of  the  6 
store  displacement  components,  enter  zero  for 
these  components. 

Format  --  (615').  Number  of  cards  is  1. 

Data  are  entered  bv  subroutine  RFADY. 

.5  i  FCALEW  Factor  to  be  used  to  scale  store  weight  prior 

to  search,  lb. 

SCAI.EI (K) ,K=1 ,3  Factors  to  be  used  to  scale  store  moments  of 

inertia  about  x,  y,  and  z  axes,  lb.-in.“ 

SCALER (K) ,K=1 , 3  Factors  to  be  used  to  scale  store  center- 

of-gruvity  x,  y,  and  z  offsets  relative  to 
nominal  store/pvlon  attachment  point,  in. 

SCAI.F.F  (K)  ,  K=  1 , 6  lectors  to  be  used  to  scale  flexibility-matrix 

diagonal  elements  for  translations  in  x,  y, 
and  z.  axes  at  nominal  storc/pylon  attachment 
point,  in.  /lb.  or  rad.  /In. -lb. 

A  suggested  approach  hr  selecting  the  above  factors  is  to  set 
each  of  them  equal  to  the  difference  between  the  maximum  and 
minimum  values  (of  the  err  respond  Ing  store  parameter)  to  be 
considered  during  the  search.  Tin’s  will  result  .in  the  store¬ 
s'*  arch  region  being  scaled  to  fit  within  a  unit  mil t i -dimensional 
-al t,  which  1 x  desirable  for  efficient  searches. 

p  trainer. er  nav  be  held  fixed  during  a  search  b*  setting  the 
ft.- responding  scale  i  n  tor  equal  to  zero. 

rermat  =  ( F.  15.5/  f  F 1  5  .  5  /  3E I  5 . 6  ,'6F  !  0 . 10  .  Number  of  cards  is  4. 


Data  are  entered  bv  ubri  et  inn  HEADY. 


■:  «  ».■  w  p>  v  m  ■ 


'A 'A  ’A  A"  a  'MVA'  T  «  V."  J"  /.»  J’ 


« 

r> 

K  • 


c 


F 


►  ■  ■ 
k .  ■ 
h  ■. 


1  r  els  8 J 


:ti  (.-k. 


ISTEP 


KmNST  =  0 


1 


Search  step  number  at  which  present  submittal 
begins . 

For  starting  a  new  search,  let  TSTEI’  =  1. 

For  a  restart  from  a  previous  search,  set 
ISTEP  equal  to  the  cycle  number  obtained  from 
the  end  of  the  printed  output  for  the  previous 
series  of  search  stpps. 

If  ISTEP  >  1,  a  separate  restart  data  file, 
obtained  via  unit  48  in  the  previous  sub¬ 
mittal,  must  be  provided  as  input  to  this 
submittal  via  unit  47,  See  Section  4. 

Store  pitch  and  yaw  inertias  are  slaved 
together  and  considered  as  a  single  search 
var iahl e . 

Using  this  option  requires  that  each  pair  of 
initial  values  and  scale  factors  for  these  two 
quantities  be  identical. 

Store  pitch  and  yaw  inertias  are  considered  as 
independent  search  variables. 


Format  =  (215, 5X).  Number  of  cards  is  1. 
Data  are  entered  by  subroutine  READY. 


VNEW 


VS 


For  a  restart  submittal,  set  VNEW  to  flutter 
speed  obtained  from  next-to-last  search  step 
in  previous  submittal  in  which  a  reduction  in 
flutter  speed  was  achieved,  knots  eas. 

For  starting  a  new  search,  set  VNEW  to  any 
large  value  that  is  greater  than  VS  below. 

Factor  to  be  used  to  scale  flutter  speed, 
knots  eas. 

Suggest  a  value  approximately  equal  to  the 
largest  anticipated  flutter  speed  for  any 
store  configuration. 


Format  =  ( E 1 0 . 1 , l OX , E I  0 .  ) ) .  Number  of  cards  is 
Data  ire  entered  by  subroutine  RF.ADY. 

Blank  card 


-4 


5-1 


1  t  em  8M 


IT0C1  =  0  Boundaries  of  store-parameter /pylon- 

flexibility  space  to  be  searched  are  defined 
by  discrete  points.  This  option  is  applicable 
only  to  tvpes  of  search  spaces  defined  in  Item 
811. 

=  1  Search-space  boundaries  are  defined  by  linear 

constraint  equations.  See  Eq .  (3-1),  page 
3-1 ,  Reference  4  . 

Format  =  (15).  Number  of  cards  is  1. 

Data  are  entered  by  subroutine  INCONS. 

Ire;1;  Logic  Item  -  No  Data 

If  IT0C1  =  0,  go  to  Item  8T.  If  ITOC1  =  1,  continue  with  items  8F 

-  8S  below. 

Item  81‘  MSTAR  Number  of  constraint  planes  (linear  constraint 

equations)  used  to  def ine  search-space 
boundaries.  Maximum  number  is  50. 

Format  --  (15).  Number  of  cards  is  1. 

Data  are  entered  by  subroutine  INCONS. 

: tem  80  Logic  Item  -  No  Data 

Repeat  items  8R  -  8S  for  I  =  1, MSTAR. ^ 

•  -  err  8K  G(.i.I)  ,J=1  ,NPARM  Each  column,  I,  consists  of  components  of  a 

unit  normal  vector  pointing  outward  from  the 
scaled  search  space  to  the  T.th  constraint 
plane.  In  general,  rows,  J,  should  correspond 
in  order  and  total  number  (NPARM)  to  the 
search  variables  for  which  non-zero  scale 
factors  are  entered  in  Item  81;  however,  if 
KCONST  =  0,  only  one  element  should  be  entered 
for  each  pair  of  pitch  and  yaw  inertias  to  be 
searched.  The  order  in  which  columns,  I,  are 
associated  with  constraint  planes  is 
arbitrary.  Maximum  total  number  of  search 
variables,  NPARM,  is  35. 

Format  -  ( 3 (E l 5 . 5 , 5X) ) .  Number  of  cards  is  (NPARM-l)/3  +  1. 

Data  are  entered  bv  subroutine  INCONS. 


'V1.;  :  r  i  nal  discussion  of  procedures  for  defining  constraint  planes,  and , 
•t-  -  r  spec  if  ically ,  Items  8R  and  88,  is  contained  in  Subsection  5.3  and 
■- ,  .  -kI  i  x  C. 


5-13 


B  ( I ) 


Scaled  distance  from  origin  to  Ith  constraint 
plane  in  direction  of  unit  norma]  vector. 
Constitutes  Ith  element  in  h  vectoi  in  Eq . 
(3-1),  page  3-1,  Reference  4. 

Format  =  (F.15.5).  Number  of  cards  is  1. 

Data  are  entered  by  subroutine  1NC0NS. 

Logic  Item  -  No  Data 

Tf  TT0C1  =  1 ,  go  to  Item  9 .  If  1  TOC  1  =  0  repeat  Items  8U  -  800  - 
for  each  store  station.  Number  of  store  stations  is  defined  by 
the  last  value  of  1DSTR  In  Item  8K. 

ITOC  =  1  Search  variables  for  this  store  station 

consist  of  store  mass,  store  pitch  or  pitch/ 
yaw  inertia,  and  up  to  six  pylon-flexibility 
parameters.  This  option  is  intended  primarily 
for  store  stations  with  single-store  ejector 
tacks . 

=  2  Search  variables  for  this  store  station 

consist  of  store  mass,  store  pitch  or  pitch/ 
yaw  inertia,  store  center-of-gravi ty  x  offset, 
and  up  to  six  pylon-flexibility  parameters. 
This  option  is  intended  primarily  for  store 
stations  with  multiple-store  ejector  racks. 

=  3  Search  variables  for  this  store  station 

consist  exclusively  of  up  to  six  pylon- 
flexibility  parameters. 

Format  =  (15).  Number  of  cards  is  1. 

Data  are  entered  bv  subroutine  1NC0NS. 

Logic  Item  -  No  Data 

If  I  ToC  =  l  and  at  least  one  of  the  first  four  factors  in  Item  81 
is  nonzero  for  the  store  station  currently  being  considered, 
continue  with  items  R.W  -  87.  below.  ii  ITOC  =  1  and  all  of  tbe.se 
factors  are  zero,  continue  with  Items  8F.E  -  8G0 . 

If  ITOC  =  2  and  at  least  er.e  of  the  first  seven  factors  in  Item  81 
is  nonzero  for  the  store  station  currently  being  considered, 
continue  with  Items  BAA  -  800.  If  ITOC  =  2  and  all  of  these 
factors  are  zero,  continue  with  Items  8EE  -  8GG. 

if  ITOC  =  3,  continue  with  items  8EF.  -  EGG. 


eir  8W  NCORN  Number  of  vertices  defining  polygon  for 

two-dimensional  store-parameter  search  space. 
Maximum  number  is  10. 

Format  =  (15).  Number  of  cards  is  I. 

Data  are  entered  by  subroutine  1NC0NS. 

t  rn  8N  logic  Item  -  No  Data 

Repeat  Item  8Y  below  for  each  of  the  NCORN  vertices.  The  vertices 
are  defined  sequentially  In  a  clockwise  direction  assuming  store 
mass  Is  the  abscissa  and  store  Inertia  is  the  ordinate  of  the 
two-dimensional  store-parameter  space. 

- i  Rv  CORNM(J)  Store  mass  at  Jth  vertex,  lb. 

CORN r (J)  Store  p^tch  moment  of  inertia  at  Jth  vertex, 

lb . -in . “ 

Format  =  (2 (K 15 . 5 , 5X)  )  .  Number  of  cards  is  I. 

Data  are  entered  by  subroutine  TNCONS. 

>r  Logic  Item  -  No  Data 

continue  with  items  8E1£  -  8GG. 

8 -I.'.  NC'ONST  Number  of  constraint  planes  defining  three- 

dimensional  store-parameter  search  space. 
Maximum  number  determined  by  requirement  that 
total  number  of  constraint  planes  be  less  than 
or  equal  to  50  considering  all  store  stations 
and  both  store-pat ameter  and  pylon-flexibility 
.search  variables. 

Format  -  (15).  Number  of  cards  is  1. 

Data  r re  entered  by  subroutine  TNCONS. 

i  .:'!■!  1 «;  ;t?rn  —  \’o  Dnt. a 

u t  r£enr?;  HC>'.  -  8DD  !'p!ow  for  each  of  the  %’•  0\ST  tr  1  :nt 


I 


:  8CC 


•  ->  m  b  D  D 


Logic  Item  -  No  Data 

Repeat  Item  8DD  below  for  three  points  on  each  constraint  plane. 
The  points  should  be  entered  sequentially  In  a  clockwise  direction 
Looking  in  toward  the  region  to  be  searched,  and  assuming  that 
store  mass,  store  pitch  inertia,  and  store  center-of-gravlty  x 
offset  (in  that  order)  constitute  the  x,  y,  and  z  axes, 
respectively,  of  a  right-handed  coordinate  system. 


CORNM(K) 

CORNI(K) 

CORNS (K) 


Store  mass  at  Kth  point  on  constraint  surface, 
lb. 

Store  pitch  moment  of  inertia  at  Kth  point  on 
constraint  surface,  lb. -in. 

Store  center-of-gravi ty  x  offset  at  Kth  point 
on  constraint  surface,  in.  Positive  for  store 
cencer-of-gravity  location  forward  of  pylon 
attachment  point. 


Format  =  (3 (El S . 5 , 5X ) 1 .  Number  of  cards  is  1. 
Data  entered  by  subroutine  INCONS. 

NSTAK 


Number  of  constraint  planes  for  pylon-flexi- 
Lil.ity  search  variables.  Only  an  upper  and  a 
lower  limit  for  each  pylon-flexibility 
variable  may  be  specified.  Maximum  number  is 


Format  =  (15).  Number  of  cards  is  1. 

Pat  3  are  entered  by  subroutine  INCONS. 
logic  Item  -  No  Data 

Repeat  Ttem  8(,G  below  for  each  of  the  NSTAR  pylon-flexibility 
constraint  planes.  It  NSTAR  -  0,  skip  Ttem  8GG. 


[  Oi  ‘K 


i  *17 


Felon-flexibility  dynamic  degree  of  freedom  to 
h"  constrained.  Enter  with  either  a  positive 
or  negative  sign  to  denote  whether  associated 
constraint  value  (next  parameter)  is  an  upper 
or  lower  flexibility  limit. 

Maximum  or  minimum  flexibility  tor  degree  of 
freedom  specified  in  previous  parameter, 
in. /lb.  or  rad. /in. -lb. 


Format  -  (T5,F)5.5).  Number  of  cards  is  i 
L.-ta  are  entered  hv  subroutine  INCONS. 


,-lb 


terns  y , 10  FASTOP  logic/data  items  not  used  in  ESP. 


t  t'.ms 

r  ems 
S 

t:  m, 


11  FASTOP  data  item  not  used  in  ESP.  (Superseded  by  separate  file 
containing  dynamic  mass  matrix.  See  Subsection  5.1) 

12-  FASTOP  logic/data  items  not  used  in  F.SP. 

16  Logic  Item  -  No  Data 

If  a  tree-free  wing  is  being  analyzed  (KH'F(37)  =  37),  continue 
below  with  Items  16A  -  17.  Otherwise  (KH!L(37)  =  0),  omit  these 
items . 

|e<\  NPMX'F  Number  of  plug  dynamic  degrees  of  freedom. 

Should  correspond  to  the  number  of  modes  in 
the  rigid-body  displacement  matrix,  and  the 
number  of  rows  and  columns  in  the  plug  mass 
matrix,  both  of  which  are  entered  In  separate 
files.  See  Subsection  b.i.  Maximum  value  is 


Format  =  (15).  Number  of  cards  is  1. 

Data  are  entered  by  subroutine  FFKASS . 

FASTOP  data  item  not  used  in  ESP.  (Superseded  by  separate  file 
containing  plug  mass  matrix.  See  Subsection  5.1.) 

FASTOP  logic/data  items  not  used  in  ESP. 


N HOOTS  Number  of  normal  modes  of  vibration  to  be 

computed.  Maximum  number  Is  40  if  Ki.lU'tiB)  * 
0.  and  (AO  -  NPGDOF)  if  KI.0E(38)  =-  it. 

NDOFFr  Number  of  displacements  per  mode  which  are  to 

he  saved  for  use  in  the  flutter-analysis 
module.  Note  that  only  out-of-plane  dis¬ 
placements  are  used  in  A  FAN  except  when  the 
elastic-axis  option  is  specified.  Maximum 
number  is  NDYDOF. 

do  Number  of  additional  displacements  per  node 

which  are  to  be  set  to  zero  for  use  in  the 
flutter-analysis  module.  See  it*n  .’1  below. 
Maximum  number  is  (220  -  NDOFF I ) . 

r:r>  :  Number  of  cards  is  1. 

of  o  r  e  enfetor!  hv  subroutine  EIGEN. 


Lwj;l.c  1 1  fin  -  No  Data 

If  ligid-body  ramies  are  not  t<-  be  included  in  the  output  irom  the 
vihration-analvsls  module  to  the  flutter-analysis  module  (KLUF(38) 
=  0),  omit  the  following  item. 

KFQ(N) ,N=1 .NPGDOF  Zcrt -airspeed  frequencies  to  be  assigned  to 

rigid-body  modes  in  flutter  analysis,  Hz.  To 
avoid  numerical  problems,  values  above  0.005 
Hz.  are  suggested;  values  less  than  0.002  Hz. 
will  be  reset  to  0.002  Hz.  within  the  program. 
See  Section  2  for  additional  comments  on  the 
use  ol  this  Input -data  Item. 

Formal  =  (3E1P.3).  Number  ol  cards  is  I  . 

Data  are  entered  b"  subroutine  V1BTFO. 

logic  item  -  No  Oat a 

The  two  variables  in  the  following  item  permit  the  user  to  reorder 
and/or  eliminate  the  modal  <i  i sp 1 acements  calculated  in  the 
vibration-analysis  module  in  preparation  for  their  use  in  the 
flutter-analysis  module.  Further,  additional  modal  displacements 
may  be  set  to  zero.  The  requirements  for  ordering  modal 
dl sp’ acements  for  the  flutter-analysis  module  are  given  in 
Reference  5,  Volume  FT,  pages  207  and  208. 

The  two  items  should  he  entered  in  pairs  1  or  1--1  .NPOFFF. 

IDFVfl'.i  Dev  ree  -of -f  reedom  number  used  in  the  vibration 

analysis  which  corresponds  to  the  flutter- 
analysis  degree-of-f reedom  number  entered  for 
the  following  variable.  Maximum  value  is 
MDYDOF. 

il)FF(I)  IVg  re  e-oi  -  r  reedom  number  to  be  used  in  the 

f  liitu'i  analysis.  Must  be  specified  in 
ascending  order.  Omitting  a  degree-of-f reedom 
number  will  result  in  the  modal  displacements 
for  that  degree  of  freedom  being  set  to  zero. 
Maximum  value  is  (NPOFFF  4  NZKRO)  . 

i  su  -  (5(214)).  Number  of  curds  is  ( (NDOFF,;- 1 ) /5  +  1). 

I'.’tn  arc  entered  by  nub  re.-.:  t  jm-  VIB1F0. 
identical  t->  FASTOP. 


'-18 


Data  Fntered  Via  Flutter-Anaivsi s  Module 


(See  also  Reference  1,  Volume  TT,  pages  .'’1  7-274,  or  Reference  3, 
Volume  II,  pages  245-302) 


Identical  to  FASTOP. 


I.C(l)  =  -1 
=  0 
=  1 


L  C  C  2 ) 


i.o(3) 


Store  search  or  p-k  flutter  analysis. 

Pressure  calculations  only.' 
k  flutter  analysis 
Divergence  analysis.' 

Number  of  vibration  inodes,  including  rigid- 
body  modes,  to  be  used  in  flutter  analysis. 
Maximum  numher  is  40. 

Number  of  lifting  surfaces.  For  the  doublet  - 
lattice  method,  which  is  the  only  aerodynamic 
theory  that  is  operational  in  ESP,  the  maximum 
numher  is  30. 


2  ( 4  ) 


h)  -=  r, 


Number  of  reduced  velocities  for  which 
aerodynamic  pressures  and/or  generalized 
forces  are  to  be  computed  or  interpolated 
depending  on  the  value  oi  EC (l 3).  If  l.C  ( 1 )  = 
0  or  l,  LCia)  must  be  less  than  or  equal  to 
30.  Tf  LC  ( 1 )  =  2,  let  LC(4)  =  1.  If  l.C  (1 )  = 
-l,  LC/4)  is  the  number  of  reference  reduced 
velocities  (see  Item  19A),  and  must  be  less 
than  or  equal  to  15. 

Number  of  air  densities  for  which  t^e  flutter 
or  divergence  analyses  will  be  run , ‘  Maximum 
numher  is  10.  If  l.C ( 1)  =  0,  let  LC(5)  =  0. 

If  KLh’F.(')  =  7,  let  LC(5)  =  1. 

Dummy  variable  in  F.SP, 

hist  calculated  pressures.' 

No  displav. 

I .  i  s  !  lilt  and  mi 'iron  t  i  ■  !  i  >  i  >  -n  t  ••  . 


I  ■  p 1  ay  . 


•-  o'  erxed  in  FSP 


LC (9)  =  1 


Krequencv-independent  additions  to  the  aero¬ 
dynamic  matrix  Q  are  to  be  read  as  data. 


b 


OUTPUT  DESCRIPTION 


Ince  the  current  version  of  ESP  was  developed  as  a  pilot  code  to  inves- 
und  demonstrate  the  store-search  procedure  in  an  expedient  manner,  the 
I  rum  the  new  code  is  not  always  as  well  structured  and/or  annotated  as 
he  desirable  in  a  production  program.  To  partially  compensate  for  this 
tier,  and  also  to  provide  an  introduction  to  the  output  from  ESP, 
ns  of  typical  ESP  output  listings  are  presented  and  discussed  in  this 
n.  Output  related  to  NASTRAN  dynamics-model  input  matrices  is  discussed 
end i x  A. 

.  ;  PRATT  ON -ANALYSIS  MODI'!  E 

•  ner  than  the  optional  Listing  of  the  input  flexibility  matrix,  the 

r  ur-  of  output  from  an  ESP  search  run  consists  of  store-parameter  data 
us t rated  in  Figure  6-1 .  NUMBER,  at  the  top  of  the  figure,  identifies 
’te— station  number;  this  value  is  the  same  as  IDSTR  that  was  input  as 

The  next  line,  DYNAMIC  DOF,  identifies  the  dynamic-degree-of-f reedom 
.<t  the  nominal,  store/pvlon  attachment  point;  these  values  are  the  same 
.DYDi'K's  tu, at  were  in;  it  as  Item  8H .  The  store-parameter  values,  which 

•  ■  1 1 e  the  bulk  of  the  output  in  Figure  6-1,  are  listed  in  the  same  order 

the  same  units  as  the  input-data  quantities  described  in  Ttems  80  - 
:X.  ''EY ,  and  PKZ,  and  PFXX,  PFYY ,  and  PFZZ  denote  pylon-flexibility 
".•sponding  to  translations  in  the  x,  y,  and  z  directions  and 
rvit  the  >. ,  y,  and  ?.  axes.  For  the  first  analysis  cycle  in  each 
~ , r  r  ,  the  INITIAL  and  CURRENT  values  are  identical,  and  the  DELTA  is 


11  it.  i,-  ;ittcr  t  lie  store-parameter  data,  but  only  for  the  first 
,  *  •  i e  In  each  ioh  submission,  is  a  listing  of  the  coefficient  matrix, 
,  vector ,  h ,  In  Eq .  (')-!),  page  3-1,  Reference  A.  When  IT0C1  =  1  in 

; C  input  data,  this  listing,  illustrated  in  Figure  6-2,  is  ident- 

tnptir  data  specified  In  Items  8R  and  8S  (except  that,  for  input, 

,  . .  f  1  is  specified  t'v  r.iws ,  whereas,  for  output,  C.  itself  is 

r  s  '.  When  i  TOC  I  =  0 ,  (.  and  1>  are  calculated  by  the  program  based 

;•  ••  r  ameter  -l.it  a  values  i  ntcred  via  Ttems  81'  -  8C.C  , 
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Figure  5-2  -  Illustrative  Primary  Input-Data  File  (4  of  4) 
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Figure  5-2  -  i llustrative  Primary  Input-Data  File  (3  of  4) 
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Figure  5-2  -  Illustrative  Primary  Input-Data  File  (1  of  4) 


.  5  -  I, I  STING  OF  SAMPLF.  PRIMARY  INPUT-DATA  FILE 


An  input-data  file  illustrating  the  use  of  many  of  the  data  items 
.inscribed  above  is  shown  In  Figure  5-2.  The  case  shown  is  for  a  maximum  of 
eight  store-parameter  search  steps  in  a  six-dimensional  space  in  which  store 
weight,  longitudinal  center  of  gravity,  and  pylon  flexibility  are  allowed 
to  vary  at  one  store  station,  and  store  weight,  longitudinal  center  of 
gravity,  and  slaved  pitch/yaw  inertia  are  to  be  varied  at  a  second  station. 

I  he  constraint  planes  which  constitute  the  search-space  boundaries  are 
let  ir.ed  in  terms  of  unit  normal  vectors  and  distances  from  the  search-space 
- 1  i i  n  , 


Note  that  each  constraint -plane  unit  normal  vector  (Item  8R)  has  at  most 
v  components,  and  that  both  components,  when  two  are  present,  pertain  to 
r.ar.eters  at  the  same  store  station.  Two  significant  implications  of  these 
characteristics  of  the  constraint-plane  data  are  as  follows.  First,  the 
".rrrLs  ible  range  of  store-parameter  variations  at  one  station  is  fully 
'-.•u-pendent  ot  the  values  at  the  other  station.  Second,  each  constraint  plane 
w : 1 1  appear  as  a  line  in  at  least  one  two-dimensional  space:  For  a  unit 

■rsa!  vector  having  two  nonzero  components,  this  two-dimensional  space  is  the 
containing  the  two  nonzero  components .  When  the  unit  normal  has  only  one 
.  :  rcnooneti  ,  the  constraint  plane  becomes  a  line  in  any  two-dimensional 

snaor?  lor  which  one  of  the  variables  corresponds  to  the  nonzero  component;  by 
'  ;  i  c '  c-n .  the  axis  corresponding  to  the  single  nonzero  unit-normal  component 
I  normal  to  this  line. 

\  .’a  discussed  In  Appendix  f,  the  specification  of  constraint  planes 

linens  ior.ai  searches  in  terms  of  two-component  unit  normals  not 
■  •  - 1  *  !  v  rnr  l  I  I  fates  tlv  def  in!  t  ior  d  the  search  boundaries,  hut  also 
i  .  ’  r'  >  an :  j  c.  1  for  producing  a  con,  trained  search  spare  with  little  open 

■■■  •  v  fne  periphery.  ’therefore,  this  approach  Is  recommended  whenever 
in  i  tnree -parameter  space  at  each  store  station  is  deemed  useful. 
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for  the  next  analysis  cycle.  However,  since  it  is  in  the  search  module  that 
ESP  achieves  its  primary  technical  distinction  from  FASTOP,  the  name  of  this 
module  has  been  changed  herein  from  that  of  its  FASTOP  predecessor. 


Nevertheless,  since  some  of  the  input  associated  with  the  FASTOP  flutter 
optimization  module  is  still  present  in  the  current  pilot-code  version  of  ESP 
rhe  practice  of  relating  ESP  input  data  to  corresponding  FASTOP  input-data 
Item  numbers  will  be  retained  in  this  module  as  well. 

item  1  i.ogic  Item  -  No  Data 

If  a  search  is  to  be  performed  to  determine  store-mass  and/or 
py lon-f lexibil itv  parameters  correspond! ng  to  a  minimum  flutter 
speed  (KLUE(7)  --  7),  continue  with  Items  2-8  below.  If  only 
an  analysis  is  to  be  performed  (KLUE(7)  -  0),  omit  these  items. 

Items  2,3  Identical  to  FASTOP 

Item  4  Blank  card  in  ESP. 

Item  5  FASTOP  .logic  item  not  used  in  ESP. 

item  6  VDES  Velocity  greater  than  maximum  anticipated 

flutter  speed  during  search,  knots  equiva¬ 
lent  airspeed. 

Format  -  (F10.3.20X).  Number  of  cards  is  1. 

Data  are  entered  by  subroutine  AFOM. 

item  7  NFIX  Maximum  number  of  search  steps  to  be  taken 

in  this  job  submission.  The  program  will 
stop  in  less  than  NFIX  steps  if  a  minimum 
flutter  speed  is  located. 

Fermat  =  (3X,I3).  Number  of  cards  is  1. 

Data  are  entered  by  subroutine  AFOM. 

“en.  S  hi  an’-*  card  in  F.SP. 
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Item  29A  Format  =  (215).  Number  of  cards  is  1. 

(  cent . ) 

Data  are  entered  by  subroutine  FI, INFO. 

Items  30-51  Identical  to  FASTOP. 

Item  5  1 A  logic  Item  -  No  Data 

If  modes  having  lew  ratios  of  on-diagonal  generalized  force  to 
corresponding  genetalized  mass  are  to  bo  automat ical 1 v  excluded 
from  a  p-k  flutter  anal  vs  is  (1(1(1)  -  -1  and  !C(38)  -  1),  enter 
data  for  the  following  item.  Otherwise  (i.C(l'i  f  or  (X  ( 38)  - 
0';,  omit  this  item. 

Oration  is  advised  when  using  tills  mode-elimination  option  in 
conjunction  with  store-parameter  searches.  Derivative 
calculations,  and  resulting  search  steps,  can  be  much  more 
sensitive  than  the  flutter  charact  eristics  themselves  to  the 
system- idea  1 izat  uu.  changes  resulting  from  mode  eliminations. 

1  ■  .;v  o  ’  PD  V  Kat  i  •»  ol  modulus  ol  on-diagonal  generalized 

force  to  value  of  corresponding  generalized 
mass  below  which  a  mode  will  be  excluded  from 
p-k  flutter  analysis. 

■:<  PV  uinal  velocity  at  which  generalized  forces 

used  In  above  ratio  are  to  be  computed,  knots 
tiis.  (Generalized  forces  are  actually 
■omnut'd  and  used  at  0. 75*V()PW,  VQDW,  and 
; .  ?r.*Vpnw. ) 

'•'•■tr.u-  (?K!D.  ’•)  Number  of  cards  is  1  . 

into  are  entered  by  suin’.  >uf  inn  FMNK5. 

•  ■:  •  d  tunica  1  f  KAt/iOP. 

•  !)k  Idert  1  .•••'.!  to  FAvrop. 

;m  SI  tip  1  >  gic/dot .!  i  t  :  nm  used  ir.  ESI*. 

■  ■  .  Si  o;:  o  f  i  ’  t  gn..  in  1.S1'. 


■  r  1  r  >  V  i  Sea  r.  h  M.  .  c.t 


P 

t  1  *  i 


1  .  1  lure  I  1  ,  pa,;  a  2/3-2  7b,  or  Reierence  3, 
ill-  ai'.  1 

1  .t  .  !i  iff. a  d  1 1  I  .*  .a  !  •;  similar  to  the  flnttor- 

1  1 '■■■  ext  ■  ■  a '  *'ia:  both  modules  determine 

■  rr  t  r  -t  ;■»  :  •  1  i  nr,  !  ■  r:  ■  r.  t  e.i  to  dr  tine  the  data 


Item  20 


RVBO(  I)  ,  I«=l  ,NRVBO  Reference  reduced  velocities. 


The  values  should  be  input  in  ascending  order, 
and  should  span  the  entire  range  of  reduced 
velocities  required  for  the  flutter  analysis. 

For  the  k  method,  this  implies  that  RVBO(l)  < 
VB0(1)  (see  Item  17),  and  that  RVBO(NRVBO) 

>  VB0(LC(4) ) . 

For  the  p-k  method,  the  following  approxi¬ 
mation  is  suggested: 

-  RVBO(l)  <  1.69*1 2*VMIN/(BR*6.28*FMAX) 

-  RVBO(NRVBO)  >  1 . 69*12*VMAX/ (BR*6. 28*FMIN) 
where 

VMIN  =  VI,  knots  (see  Item  19) 

-  VMAX  =  VI  +  (NV-1 ) *DV ,  knots 

-  FMAX  and  FMIN  are  the  maximum  and  minimum 
zero-airspeed  frequencies  in  Hz. 

-  BR  is  the  reference  semichord,  in.  (see 
Item  1 5) . 

If  FMIN  <  0.0 l,  let  RVBO(NRVBO)  =  1.0E+05  and 
RVBO(NRVBO-l)  »  2000. 

Format  =  (7E10.0).  Number  of  cards  is  ((NRVB0-l)/7  +1). 

Data  are  entered  by  subroutine  FLINFO. 

Item  21-25  Identical  to  FASTOP.1 
Item  26  Logic  Item  -  No  Data 

If  no  structural  damping  Is  added  to  the  stiffness  matrix  (LC(16) 

=  0),  omit  the  following  three  items,  and  go  to  Item  29A. 

If  different  structural  damping  values  are  added  to  the  complex 
stiffness  matrix  for  each  mode  (LC(16)  =  -1),  omit  Ttem  27,  and  go 
to  Item  28. 

If  the  same  value  of  damping  is  added  for  all  modes  (LC(16)  =  1), 
enter  data  for  Item  27  and  omit  Items  28  and  29. 

Item.;  27-29  Identical  to  FASTOP. 


Item  :NA  NRBTR 


Number  of  rigid-body  translation  modes 
Included  in  IFLMD(I),  Item  14. 


NK  i’TOT 


Total  number  of  rigid-body  modes  included  in 
I FLMD ( I ) . 


Not  applicable  to  store-parameter  search  runs. 


1  cem  16 


Logic  Item  -  No  Data 


For  divergence  analysis  (LC(1)  =  2),  omit  Items  17  -  34,  and  go  to 
Item  35. 

For  steady-state  pressure  calculations  only  (I.C(l)  =  0  and  LC(33) 

-  1),  omit  Items  17  -  55,  and  go  to  Item  56. 

For  n-k  flutter  analysis  (LC(1)  =  —  1 ) ,  omit  the  following  two 
'terns,  and  go  to  item  19. 

For  k  flutter  analysis  (I.C(l)  =  1),  or  for  oscillatory  pressure 
calculations  (LC(1)  =  0  and  I,C(33)  =  0),  continue  with  Item  17 
be  low . 

id’  f  1 )  ,1  =  1  ,1.C  (  * )  Reduced  velocities  to  he  used  in  k  flutter 
analysis  or  in  oscillatory  pressure 
calculations. 

F-rnvt  ' , 'Ll  0.0) .  Number  of  cards  is  ((LC(4)-i)/7  +  1). 

i  i'o  ire  entered  hv  subroutine  FI. INFO. 

u  '  !  r-m  -  No  'Inca 

T  t -  in  In,  this  i  tem  is  ranched  only  for  k  flutter  analysis 
*  -  11  m  for  oscillatory  pressure  calculations  ( I.C ( 1 )  =  0 

..d  Ia.'.U;  =  03. 

k  *  1  ii  1 1  e  analysis  with  generalized-aarodynamic-force 
. nb'rpr ’ ni ion  C.Cil)  »  1  and  LC(13)  =  1),  omit  the  following  item, 

-■  .d  xo  to  Item  !°A. 

For  k  > !ut  fer  innivsis  with  directly  computed  generalized  forces 
K'Hi  •  1  oid  l.C  i,  1  3)  -  0)  ,  omit  Items  19  -  20,  and  go  to  Item  21. 

i  or  use  i  t nrv  pressure  calculations  (LC(1)  »  0  and  LC(33)  =  0), 
ms  19  -  55,  and  go  to  Item  56. 

■  n  t  FASTOF. 

Number  of  reference  reduced  velocities,  i.e., 
'■umber  of  reduced  velocities  and  corresponding 
directly  computed  general ized-nerodynnmie- 
force  matrices  that  will  be  used  as  a  basis 
for  generalized-force  interpolation.  Maximum 
value  Is  15. 

: r  :  •  i 9  > .  Number  rf  cards  is  1. 

.ir  ci’terid  hv  subroutine  FLL.1F0. 


Item  4 
(cent . ) 


0 


No  display. 


T  teni  I 


l  :  vr-..i 


C.p*  1 

'  ee 

Not 


I,C(38)  =  I  User  will  input  a  ratio  of  on-diagonal 

generalized  force  to  corresponding  generalized 
mass  which  will  be  used  for  automatic 
exclusion  of  vibration  modes  from  a  p-k. 
flutter  analysis. 

=  0  All  modes  initially  selected  by  user  will  be 

used  in  the  flutter  analysis. 

If  LC(1)  4  -1,  let  LC ( 38)  =  0. 


Format  =  (1015).  Number  of  cards  is  4. 


Data  are  entered  by  subroutine  AFAM. 

1  3 

IN  =  1  Modal  vibration  data  are  input  on  cards.  * 

-  2  Not  used  in  ESP. 

-■  3  Modal  data  are  obtained  from  vibration- 

analysis  module. 

Format  =  (15).  Number  of  cards  is  1. 


bat  a  are  entered  by  subroutine  POOL. 
t>—  I  ?  Identical  to  FASTOP.^ 


FAS'! OP  data  item  not  used  in  ESP. 

I  FLMD(  I )  ,  I  =  1  ,LC  (2)  Indices  of  modes  from  vibr.it  ion-analysi  s 

module  to  be  used  in  the  flutter  analysis.  1 
is  the  mode  index  in  the  flutter-analysis 
module  that  corresponds  to  1FLMD  in  the 
vibration-analysis  module. 

For  mat  ~  '  1015).  Number  of  cards  is  (  (l.C  ( 2 ) - 1 )  /  1  0  +  1). 

Drii.i  are  entered  bv  subroutine  POOL. 


Identical  to  FASTOf  . 


o,i  net  checked  In  ESP. 

cautionary  comment  in  footnote  to  Item  51A. 
applicable  to  store-parameter  search  runs. 


I  t  em  4  IX f  29)  =  1 

(cont  . ) 

-  0 

ix  no)  =  i 

=  o 


ixo  I)  =  i 


=  0 

1X02)  =  I 


Display  physical  vectors  correspond ing  to 
the  displayed  modal  eigenvectors. 

No  display. 

Display  flutter  determinant  in  k  flutter 
analysis  (see  LC(12)). 

No  displav. 

If  1X(1)  =  -1  or  0,  let  LCOO)  -  0. 

User  will  input  changes  to  the  general  ized 
masses  and  the  modal  frequencies. 

No  changes. 

User  will  input  revisions  to  the  generalized 
stiffness  matrix. 


=  0 

IX  (33)  =  1 
=  0 


!  ((14)  “  I 


=  0 

1X05)  =  o 


irrilj)  =  1 


=  0 


No  revisions. 

,  1 

Steadv-state  analysis. 

Oscillatory  analysis. 

It  LC ( 1 )  =  2,  let  IX(33)  =  1. 

i'ser  will  input  factors  to  scale  the  computed 
aerodynamic  forces. 

No  factors. 

Fixed  value  in  ESP.  1X(35)  =  1  is  associated 
with  Mach-box  aerodynamics,  which  is  not 
operational  in  ESP. 

Compute  eigenvectors  and  the  aerodynamic  force 
gradients  required  for  flutter  redesign. 

Do  not  compute. 

If  JX <  1 )  -1  ,  let  IX (36)  =  0. 

if  KTXK(7)  =  7,  eigenvectors  and  gradients  are 
always  computed. 


lot  ilnuli  1  e  t  -  1  a  t  t  ice  program,  displav  geometru 
data  .to  (U  i  a  t  ed  with  basic  douhlel  elements. 


1  i  i 


i  r  em  4 
( font  .  1 


LC (22) 


0 


=  1 

I.C(2  3)  =  1 
=  0 


1.0(2  A)  =  1 

=  0 


LC (25) 


’.(•(26) 


i..r(27) 


LC ( 28)  =  1 


=  0 


Compute  AIC  arrays  as  part  of  this  job 
submission  and  save  as  output  data  file. 

AIC  arrays  exist  as  an  input  data  file,  and  do 
not  need  to  be  recomputed. 

Display  modal  input  data.' 

No  display. 

Display  interpolated  modal  data.' 

No  display. 

Number  of  user-specified  mode-elimination 
cycles  requested  for  the  flutter  analysis. 
(Distinct  from  automatic  mode  elimination  9 
specified  via  l.C(38).)  Maximum  number  is  25." 
Tf  KLU F(7)  =  7,  let  LC(25)  =  0. 

Number  of  stiffness-variation  cycles  requested 
for.the  flutter  analysis.  Maximum  number  is 
20/  If  KLUE(7)  =  7,  let  LC(26)  =  0. 

Index  of  the  vibration  mode  whose  stiffness  is 
to  be  varied  in  the  flutter  analysis.  If 
LC (76)  =  0,  let  LC ( 2 7 )  =  0. 

Display  eigenvectors.'  If  LC ( 1 )  “  -1,  the 
eigenvectors  for  the  critical  flutter  root  in 
a  user-specified  velocity  interval  is 
displayed.  If  LC(1)  =  1,  the  eigenvectors  for 
all  roots  between  user-specified  reduced 
velocities  and  real  frequencies  are  displayed. 

No  display. 

If  LC ( 1 i  -  0  or  2,  let  LC(28)  =  0. 

If  i.C(l)  =  -1  and  LC(36)  =  1,  the  eigenvector 
corresponding  to  the  critical  flutter  root  and 
flutter  speed  is  always  displayed. 


sc  ; on  not  v  becked  In  ESC. 

’  eot  i  becked  for  value  other  than  zero. 

•  :a!  saantitles  saved  are  not  AlC’s,  but  intermediate  results  of 

q  • '  •  ■.  t  ft-  •  solution  procedure  (see  Reference  7,  pages  34-35)  . 


f  1  cm  -4 
(cont.) 


I.0(1M  =  0 


No  structural  damping  is  addec}  to  the  complex 
generalized-st ( ffness  mat  rix. 


1.0  ( 1  7 )  - 

u:( 18)  = 


i.rd'i)  = 

1,0(20)  - 

1 1  ~- 


:  r  "e,t  chocked 

i  '  ’  arm  (rut  s t ■ 
a  r  ,  the  i  n  r  odu 
■  :  iiu.ne  t  f  c a  ! 

V  i  ■  V-  !  »:  V'  IS  <  >  t 


-1  Different  damping  values  are  added  to  the 

generalized-stiffness  matrix  for  each  inode. 

I  The  same  value  of  damping  is  added  for  all 

modes . 

!  Display  the  number  of  iterations  required  to 

obtain  each  root  in  the  p-k  flutter  analysis. 

0  No  displav. 

If  LC ( 1 )  t  -1  ,  let  LC (17)  -  0. 

1  For  the  third  and  higher  velocities  in  the  p-k 

flutter  analysis,  the  initial  estimate  of  each 
root  is  obtained  hv  extrapolating  from  the 
root  values  at  the  previous  two  velocities. 

0  The  value  of  the  root  at  the  previous  velocity 

is  used  as  the  root  estimate. 

If  L.C(l)  J  -1  ,  let  LC  ( 1 8)  =  0. 

1  Order  the  roots  a|'ter  solution  by  the  p-k 

flutter  analysis. 

0  No  ordering. 

If  1.0(1  )  ^  -1  ,  let  1.0(19)  =  0. 

1  Display  the  root  iterations  in  the  p-k  flutter 

analysis  (l.C(L)  =  -1),  or  display  intermediate 
results  of  the  k  flutter  analysis  (1,0(1)  - 
1). 

0  No  d ! sp 1 av . 

Fixed  value  in  ESP.  Denotes  use  of  subsonic 
douh 1 e t- 1  a 1 1 ice  aerodynamics.  Supersonic 
Much-box  and  subsonic  assunied-pressure- 
t  unction  methods  are  not  operational  in  F.SP. 


in  HSP . 

re-parameter  searches,  this  option  is  not  recommended . 
ct  ion  of  realistic  structural-dumping  values  is  suggested 
problems  that  might  he  encountered  when  some  modes  have 
<ie  rod  vnam  i  c  damp  i  ng  . 


1-2  1 


NUHBER= 


1 


DYNAMIC  DOF  =  173  174  175  176  177  178 


INITIAL  VALUE  CURRENT  VALUE  DELTA  VALUE 


WEIGHT . . 

. 28450E+04 

.  30066E  +  04 

.  16156E  +  0  3 

IXX  .  . 

. 10000E+01 

. 10000E+01 

0 . 

I  Y  Y  .  . 

. 63700E+0? 

. 63700E+07 

0  . 

L  Z  7.  .  . 

. 63700E1 07 

.  63700E  +  07 

0  . 

RX  .  . 

. 20837E+02 

.  21376E  +  02 

.  53909E  +  00 

P  Y 

0  , 

0. 

0  . 

RZ  .  . 

. 13000E  102 

.  1 3000E  +  02 

0  . 

PFX  .  , 

. 1 1 162E- 04 

.111 62E-04 

0  . 

PFY  .  - 

.39162E  04 

. 39162E-04 

0  . 

PFZ 

.53 5 53 E -04 

.  53553E-04 

0  . 

PFXX  .  . 

.  13290E--0<-. 

. 1 3290E-06 

0  . 

PFYY  .  , 

. 28086E -07 

. 30183E-07 

.  20970E-08 

PFZZ  .  . 

.  6 1 7  3  0  E  -  0  7 

. 6 1 730E-07 

0  . 

NUMBER—  2 

DYNAMIC  DOE  =  190  191  192  193  194  195 

INITIAL  VALUE  CURRENT  VALUE  DELTA  VALUE 


WEIGHT  . 

. 35000E+04 

.  35000E  '•04 

0  . 

T  X  X  .  . 

•  10000 ErO  1 

. 1 OOOOE+Ol 

0  . 

IVY  .  . 

. 95000EF07 

.  95000E  +  07 

0  . 

IZ2  .  . 

. 9 5000E+07 

.  95000E  +  07 

0  . 

R  A  .  . 

•  .  A  8 1  9  0  E  -!•  0  1 

-  .  64486E  +  01 

, 

R  Y  ,  . 

0  . 

0  . 

0  . 

R  7  ■  . 

. 1 3000E+02 

.  1  3000E  +  02 

0  . 

P  F  X  .  . 

. 47387E-05 

.  47387E-05 

0  . 

P  F  Y  .  . 

.  5  9  4  2  6  E  ~  0  4 

.  39426E-04 

0  . 

P  f  7.  .  . 

.  1  7  2  5  7  £  -  0  4 

.  17257  E. -  04 

rs 

\S  » 

tf  XX  .  , 

.  12662006 

.  12662E-06 

o  t 

P  F  .. 

.  2  7  6  9  3  £  -•  0  7 

.  27693E-07 

0  . 

PF  7.2  .  . 

.  603922-07 

.  60892E-07 

0  . 

‘yre  h  -1  -  rvplu-i.:.  Listing  of  St.  ro  -Parameter  Data  from  Vibration- 


Analysis  Module. 
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C  Matrix  -  Components  of  Unit  Normal  Vectors  Pointing  Outward  from  Search  Space  to  Each 
Constraint  Plane  (Row  Number  Defines  Component  Direction;  Column  Number  Defines  Constraint 
Plane) 


Prior  to  entering  the  eigenvalue  analysis  for  the  first  analysis  cycle, 
the  dynamics-model  mass  matrix  that  was  entered  as  a  separate  data  file  (see 
Subsection  5.1)  is  listed.^  This  output,  partially  illustrated  in  Figure  6-3, 
would  be  the  final  mass  matrix  used  for  the  eigenvalue  analysis  in  the  case  of 
an  analysis-only  run  (KLUE(7)  =  0  in  input-data  Item  6  for  the  main  program). 
However,  for  a  store-search  run  (KLUE(7)  =  7),  the  elements  in  this  mass 
matrix  that  correspond  to  the  store  degrees  of  freedom  are  replaced  subsequent 
to  the  listing  by  new  elements  calculated  from  the  data  in  Figure  6-1  accord¬ 
ing  to  their  analytical  definitions  contained  in  Reference  4,  Eq .  (4-2),  page 

2 

4-3.“  The  replacement  elements,  illustrated  in  Figure  6-4,  are  listed  immedi¬ 
ately  after  the  original  total  mass  matrix  for  the  first  analysis  cycle  in  a 
search,  and  updated  values  for  these  elements  are  listed  in  each  subsequent 
analysis  cycle  as  the  search  progresses. 

At  the  conclusion  of  the  eigenvalue  analysis,  the  frequencies  of  the 
flexible  modes  are  listed  adjacent  to  modal  indices  that  begin  with  1  for  the 
first  flexible  mode  (see  Figure  6-5(a)).  In  a  subsequent  listing  of  mode 
shapes,  rigid -body  modes  are  given  along  with  the  flexible  modes  If  KLUE(38) 
in  the  main-program  input  data  is  38  (see  Figure  6-5(b)).  Here,  the  modal- 
index  values  begin  with  1  for  the  first  rigid-body  mode,  and  the  previous 
t lexib 1  e-mode  indices  have  been  incremented  by  the  number  of  rigid-body  modes. 
This  combined-mode  indexing  system  i3  followed  throughout  the  remainder  of  the 
output  from  the  vibration-analysis  module. 

o.2  -  FLUTTER-ANALYSIS  MODULE 


A  typical  listing  from  the  first  page  of  output  from  the  flutter- 
nr.aiysis  module  is  shown  In  Figure  6-6.  Here,  the  modal  indices  are  changed 
igain,  this  time  based  on  the  modes  remaining  following  the  selection  process 
iefired  by  Item  i  cf  the  f  1  utter-ana  lys  1  s-module  input  data. 


rbis  listing  is  obtained  from  search  runs  only.  However,  an  equivalent 
listing,  produced  immediately  after  reading  the  mass  matrix,  is  provided 
bv  setting  KLL'F.V(8)  =  8. 

Modi  (feat  ’mis  to  Ref**  rent  e  h  including  Eq  .  (4-2)  ,  are  discussed  in 
Subset 1 1  on  6 .  ' . 


n  ' 
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L ; 

MATRIX  NAHE=  REF  HASS  (  206  X 
( DYNAMIC  HASS  HATRIX  ENTERING 

206)  PRINT  LOWER  TRIANGLE 

CURRENT  FOP  RUN  -  DOES  NOT  REFLECT  STORE  REDESIGN) 

• 

n. 

ROW 

COL 

VALUE 

ROW 

COL 

VALUE 

ROW 

COL 

VALUE 

;.i 

1 

l 

2.699999E+01 

* 

’  .  'J 

2 

2 

2 . 699999E+01 

..  -• 

3 

3 

2 . 699999EI01 

.  \ 

4 

4 

4 . 749995E+02 

-  - 

5 

3 

4.047360E+00 

5 

5 

2 . 060606E+03 

v/ 

w 

6 

2 

4.047360E+00 

6 

6 

2.460606E+03 

• 

] 

17  o 

168 

-1 . 923599E+03 

170 

170 

• 

2 . 952039E  i  04 

* 

< 

171 

167 

1 . 923599E+03 

171 

169 

3.984463E+02 

171 

171 

6.649581E+04 

« 

172 

1 68 

3 . 984463E+02 

172 

170 

-5. 578250E+03 

172 

172 

3 . 747545E+04 

173 

173 

9 . 999996E-02 

174 

174 

9.999996E-02 

* 

175 

175 

9.999996E-02 

J 

m 

176 

174 

-1 .519999E+00 

176 

176 

2 . 410397E+01 

177 

173 

1 . 519999E+00 

177 

175 

4.199950E-01 

177 

177 

2.586791E+01 

;  ■ 

173 

174 

4.199950E-01 

178 

176 

-6.383921E+00 

178 
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2.763956E+00 

1 

179 

179 

9 . 999996E-02 
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190 

i  80 

9.999996E-02 
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c  P 

134 

2 . 033520E+03 

188 

186 
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3.984800E+02 

188 
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6.966169E+04 

< 

g 

r. 

j.39 

165 

3 . 984800E+02 

189 

187 

-5 . 897504E+03 

189 

189 

3 . 747565E+04 

« 

1/0 

190 

9 .999996E-02 

i  9  i 

171 

9 . 999996E-02 

4 

.  ** 

1  V  \ 

i9: 

*  7 1 

9.999996E-02 

-8 . 300000E-01 

193 

193 

7 , 889000E+00 

- 

;  ?  ii 

i  •?  0 

8.300000E-01 

194 
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194 
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9 . 184927E+00 

:  .  J 

■ 
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■  s  1 
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195 
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9 . 999996E-02 
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r**  - 

l  -  Tvplc.il  Listing 

from 

m 

Vibration-Analysis  Module  of 

Dynamics- 

j 

* 

Model  Mass  Matrix  as 

Entered  via  Separate 

Data  File. 

•s 

NEW  HASS  MfiTkIX  FOR  STORE  HUMBER 


fREC.-lNClES,  ;  3  KQf.flALiZED  EIGENVECIQha  eGR  ALL  VIBRATION -ANALYSIS  DEGREES  Of  FREEDOM  (ABSOLUTE  MOTION) 
RlGiD-BOLT  MODES  (IF  REQUESTED)  FOLLOWED  BY  FLEXIBLE  HQDE:i 
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Figure  6  5  Typical  Listing  from  Vibration-Analysis  Module  of  Flexible-Mode  Frequencies  and  Rigid-Body  and 
Flexible  Mode  Shapes. 
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Figure  6-6  -  Typical  Listing  of  Modal  Data  from  Flutter-Analysis  Module. 


If  automatic  mode  elimination  based  on  ratios  of  generalized  forces  to 
generalized  masses  has  been  selected  (LC(38)  =  1  in  Item  4  of  flutter- 
analysis-raodule  input  data),  the  output  illustrated  in  Figure  6-7  is  obtained. 
The  data  provided  in  Item  5 IB,  along  with  the  two  additional  velocities 
calculated  by  the  program,  are  listed  at  the  top  of  the  figure,  followed  by 
the  results  of  the  mode-elimination  test. 

Typical  CalComp  plots  of  flutter-analysis  results  are  shown  in  Figures 
6-8  through  6-10.  The  first  of  these  figures  is  representative  of  output  from 
a  p-k  flutter  solution  that  is  done  either  during  the  first  analysis  cycle  of 
a  store-search  run  or  during  an  analysis-only  run.  The  second  plot  illustrates 
the  more  abbreviated  p-k  output  obtained  from  a  second  or  subsequent  analysis 
cycle  of  a  store-search  run.  A  k-method  flutter  solution  is  shown  in  the 
third  of  this  series  of  figures. 

6.3  -  SEARCH  MODULE 

The  output  from  the  store-parameter  search  module  is  discussed  almost  in 
its  entirety,  since  this  ESP  output  is  completely  new  compared  to  FASTOP. 

After  repeating  the  flutter  speed  (in  knots  equivalent  airspeed)  that  was 
previously  listed  In  the  output  from  the  flutter-analysis  module,  the  deriva¬ 
tives  of  the  flutter  speed  with  respect  to  the  various  store  parameters  are 
listed,  together  with  several  intermediate  results  (see  Figure  6-11).  This 
output  is  grouped  according  to  store-station  number,  with  the  derivatives 
being  listed  first  using  the  same  store-parameter  notation  that  was  discussed 
previously  in  connection  with  the  vibration-analysis  module  (cf.  Figure  6-1). 
The  three  arrays  Immediately  thereafter  are  the  matrices  in  Eq.  (4-5),  page 

4-4,  of  Reference  4;  SI,  S2,  and  S3  in  the  listing  correspond  to  S  ,  S  ,  and 

x  y 

in  the  reference.  The  U  and  V  vectors  given  next  are  the  flutter  eigen¬ 
vector  and  its  associated  row  vector  in  physical  coordinates  for  the  store 


The  title  page  for  this  module,  however,  retains  the  previous  FLUTTER 
OPTIMIZATION  MODULE  wording. 


MODAL  ELIMINATION  BASED  ON  RATIO  OF  GENERALIZED  FORCES  TO  GENERALIZED  HASSES 
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Figure  6-7  -  Typical  Listing  oi  Automatic-Mode-Elimination  Results  from  Flutter-Analysis  Module. 
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Figure  b-8  -  Typical  CalComp  Plot  of  p-k  Flutter-Analysis  Results  from 
Either  an  Analysis-Only  Run  or  the  First  Analysis  Cycle  of 
Store-Search  Run. 
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Figure  6-11  -  Typical  Listing  from  Search  Module  of  Flutter  Speed  and  Flutter-Speed  Derivatives,  Including 
Intermediate  Results  (2  of  2). 


.--on  saved  from  the  previous  cycle.  The  data  set  that  should  be  used  for  the 
i estate  is  the  last  one  that  contains  information  corresponding  to  an  analysis 
v .  i .  •  in  which  the  flutter  speed  was  decreased.  Thus,  if  the  next-to-the-last 
ina lysis  cycle  exhibited  a  lower  flutter  speed  than  the  one  before,  the  last 
j.:t  i  set  in  the  TAPE48  file  should  be  used  as  an  input  file  (via  TAPE47)  in 
the  next  job  submission. 

'.'he  TAPE48  file  is  similar  to  the  TAPE40  file  to  the  extent  that,  for  the 
a >st  part,  it  too  consists  of  selected  portions  of  the  printed  output, 
i.-wocor,  as  seen  in  Figure  6-18,  the  TAPE48  file  is  not  annotated,  and  there¬ 
to  a  brief  discussion  follows  in  which  the  lines  in  Figure  6-18  are  related 
■  T  responding  lines  in  Figures  6-12,  6-13,  6-15,  and  6-16.  All  these 
res  contain  output  from  the  same  job  and  analysis  cycle. 


;.e  first  two  lines  in  Figure  6-18  correspond  to  the  PREVIOUS  SCALED 
<h!.S  in  Figure  6-13,  and  immediately  thereafter  is  the  OLD  INVERSE 

in  Figure  6-12.  The  next  array  is  a  set  of  previous  scaled  deriva- 
I'his  array  does  not  have  an  exact  counterpart  in  Figures  6-12  and 
.  but  it  would  have  been  printed  as  the  SCALED  DERIVATIVES  in  the  version 
: ere  6-13  corresponding  to  the  previous  analysis  cycle.  The  next  two 
,  with  one  number  per  line,  are  the  step  size  and  the  flutter  speed  for 
r ; *  ions  cycle.  These  would  have  been  listed  in  the  version  of  Figure 
.  r  Chat  cycle. 


r.vxt  group  of  data  in  Figure  6-18  consists  of  the  number  of  active 
.  rn.-t  tve  constraints  at  the  end  of  the  previous  cycle,  and  corresponds  to 
mt>-.  r  I  .'ii  shown  at  the  top  of  Figure  6-12.  Next,  the  previous-cycle 
.  MoMFNi  MATRIX  in  Figure  6-12  is  listed.  The  final  line  in  Figure  6-18 
f  no  .-.ur.be  i  of  new  constraints  introduced  by  the  previous  step, 

■  o  '.I.-  It’d  !  ..-t-s  of  t  h*>.se  constraints.  Tf  the  first  number  is  zero, 

•  i;. i  number  is  the  index  ot  the  nearest  inactive  constraint  in  the 

.  -  ,  i  ’he  step.  Tills  would  have  been  listed  as  .TNEW(l)  near  the  bottom 
- •  -i  o  rf  Figure  6-13  for  the  previous  analysis  cycle. 


2>8919E-01  5 . 08220E -0 1  1.00308E+00  1.00000E+00  1.00000E  +  00 
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Figure  6-18  -  Listing  of  Typical  Restart  File  from  Search  Module 


NEW  INVERSE  HESSIAN  AND 

9.98846E-01  1.30815E-02 

DIRECTION 

1.90237E-03 

-2.28088E-02 

1.52897E  02 

-1.41320E-02 

1 . 30815E-02  8 . 51723E-01 

-2.15632E-02 

2.58536E-01 

-1.73308E-01 

1 .60185E-01 

1.90237E-03  -2.15632E-02 

9.96864E-01 

3.75975E-02 

-2.52032E-02 

2.32948E-02 

-2 . 28088E-02  2.58536E-01 

3. 75975E-02 

5.49218E-01 

3.02178E-01 

-2 . 79297E-01 

1 . 52897E-02  -1.73308E-01 

-2.52032E-02 

3.02178E-01 

7 . 97437E-01 

1 .87225E-01 

-1.41320E-02  1.60185E-01 

2.32948E-02 

-2.79297E-01 

1.87225E-01 

8 . 26952E-01 

3 . 7  4836E-02  -1.07947E-01 

3 . 40077E-03 

-5 . 46144E-02 

-4.40454E-03 

-2.52630E-02 

SCALED  DERIVATIVES  ARE 
3 . 33739E-02  -6.13634E-02 

1.01752E-02 

-1.35837E-01 

5.00426E  02 

-7 . 55873E-02 

PREVIOUS  SCALED  VARIABLES  ARE 
7.63919E-01  5 . 08220E-01  1.00308E+00 

l.OOOOOE+OO 

l.OOOOOE+OO 

-1.7 484AE-0 1 

PRESENT  SCALED  VARIABLES  ARE 
8.1258SE-01  5.21368E-01  1.Q7797E+00 

l.OOOOOE+OO 

l.OOOOOE+OO 

1 .65349E-01 

NORM  OF  SCALED  GRADIENT  VECTOR  =  1.7791E-01 


NEH  SCALED  VARIABLES  ARE 

S  >  X 1 172E-01  5.20943E-01  1.07692E+00  l.OOOOOE+OO  l.OOOOOE+OO 

-1.57530E-01 


»t*«  NEW  STORE  PARAMETERS  *ttt* 


PARA 

STORE  NO.  1 

STORE  NO.  2 

STORE  NO. 

U 

3 . 00134E+03 

3.50000E+03 

IX 

1 .OOOOOE+OO 

l.OOOOOE+OO 

ir 

2  - 37000E+06 

9.50000E+06 

IZ 

6 . 37000E+06 

9 . 50000E+06 

RX 

2 . 1 3587E+01 

-6 . 14367E+00 

R  ! 

> , 

0. 

?Z 

l  .  30000E+01 

1 . 30000E4  t-i 

F  X 

1.11620E-05 

4.73867E-06 

Ft 

o .  91 623E-05 

3 . 94256E-05 

FZ 

5.35534E-C5 

1 .72570E-05 

F  XX 

1  .  32900E-07 

1 .26623E-07 

F  '( 

3.01537E-08 

2 . 76931E-08 

1*  v.  1. 

6. 1730 4E -08 

6 . 08923E-03 

CURRENT 

STEP  SIZE  = 

1.000 

CURRENT 

FLUTTER  SPEED 

=  363,974 

KNOTS 


i  •>  6-17  -  Listing  of  Typical  Summary  File  from  Search  Module  (2  of  2) 
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DERIVATIVES  FOR  STORE  NUMBER  1 


U  9.01997E-03 
IXX  1 . 76939E-05 
IYY  -2.41730E-06 
IZZ  -2 . 41730E-06 


DERIVATIVES  FOR  STORE  NUMBER  2 
U  -3.88106E-02 

IXX  -4.55362E-05 

IYY  5 • 26764E-06 

IZZ  5 . 26764E-06 


STORE  PARAMETERS  FOR  REDESIGN  CYCLE  NO.  2 
THE  FLUTTER  SPEED  HAS  BEEN  DECREASED  BY  7.1235  KNOTS. 
PREVIOUS  MOVE 


4.3A457E-02  1.31486E-02  7.48922E-02  0.  0. 

9.49678E-03 

DELTA  GRAD 


3.74929E-02  8.31172E-02  8.50673E-02  -1.21997E-01  8.17798E-02 

-6.60905E-02 


OLD  INVERSE  HESSIAN 

l.OOOOOE+OO  0. 

0. 

0.  l.OOOOOE+OO 

0. 

0.  0. 

0. 

0.  0. 

0. 

0.  0. 

0. 

0.  0. 

l.OOOOOE+OO 

HOLD*(DELTA  GRAD) 


3.74929E-02  8.31172E-02  3.50673E-02  -1.21997E-01  8.17798E-02 

-A.60905E-02 

Z=  DELTAX-HOLDG 


6. 17284E-03  -6.99686E-02  -1.01752E-02  1.21997E-01  -8.17798E-02 

7 . 55873E-02 

CK  AND  Z2  =  - . 33017E-01  .32322E-01 


0. 

0. 

l.OOOOOE+OO 

0. 

0. 

0. 


0. 

0. 

0. 

l.OOOOOE+OO 

0. 

0. 


0, 

0. 

0. 

0. 

l.OOOOOE+OO 

0. 


Figure  6-17  -  Listing  of  Typical  Summary  File  from  Search  Module  (1  of  2) 
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Nf1-'  SCALED  variables  are 

1172E-01  5.20943E-01  1  .  C7692E  +  00  l.OOOOOE+OO  l.OOOOOE+OO 

7530E-0I 
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Figure  6-16  -  Typical  Listing  of  Search-Module  Output  Including  New  Store  Parameters  in  Terms  of  Physical 
Units 


Below  PG,  the  inactive  constraint  indices  for  the  current  search  point 
are  isted,  followed  by  several  quantities  defining  the  next  search  step. 
JNEW(l)  is  the  index  of  the  nearest  inactive  constraint  plane  in  the  direction 
of  the  search  step.  STEPS  is  the  step  size  required  to  reach  that  constraint, 
and  STEPSS  is  the  minimum  of  the  two  quantities  STEPS  and  1 .  VALU  is  equal  to 
0.33/PNORM,  and  STEP,  which  is  the  final  step  size,  is  the  minimum  of  the  two 
quantities  VALU  and  STEPSS.  A  further  discussion  of  the  step-size  determina¬ 
tion  procedure  is  given  in  Reference  4,  Subsection  3.4,  page  3-9  and  3-10. 

In  Figure  6-16,  the  new  scaled  variables  have  been  converted  to  physical 
store-parameter  values.  Also,  the  scaled  step  size  and  the  value  of  the 
flutter  speed  before  the  new  step  are  repeated. 

This  concludes  the  discussion  of  the  printed  output  usually  obtained  from 
the  search  module.  However,  as  implied  in  this  discussion,  additional  output 
will  be  printed  in  some  situations,  such  as  search  steps  that  require  a  line 
search  (see  page  3-10  of  Reference  4),  that  involve  a  move  off  a  constraint 
plane,  or  that  terminate  at  a  local  minimum. 

In  addition  to  the  printed  output,  the  search  module  also  provides  two 
formatted  disk  files  via  TAPE40  and  TAPE48  as  discussed  in  Section  4.  The 
first  file  contains  selected  portions  of  the  printed  output,  including  annota¬ 
tions,  and  is  intended  to  facilitate  a  quick  review  of  progress  made  during  a 
search,  in  preparation  for  a  restart  of  the  search  if  necessary.  A  comparison 
of  the  illustrative  example  in  Figure  6-17  with  the  output  in  the  previous  six 
figures  provides  a  complete  description  of  its  contents.  Although  not  shown 
In  the  figure,  the  summary  material  in  a  TAPE40  file  includes  results  for  all 
the  analysis  cycles  executed  in  the  run  from  which  it  is  generated. 

The  file  obtained  via  TAPE48  consists  of  a  series  of  results,  such  as 
those  shown  in  Figure  6-18,  that  car.  be  used  to  restart  a  search  in  a  manner 
that  takes  advantage  of  the  previously  computer  inverse-Hessian  matrix.  A  new 
'■c  of  potential  restart  data  is  obtained  from  each  analysis  cycle  except  the 
1 irst,  and  the  information  that  is  written  in  each  cycle  is  data  that  has 
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Figure  6-15  -  Typical  Listing  of  Search-Module  Scaled  Output  Defining  Constrained  Search  Step, 

Convergence-Test  Parameters,  and  Step-Size  Parameters. 


inverse-moment  matrix  is  calculated  and  listed  once  to  reflect  the  new  inverse 
Hessian  matrix;  this  is  the  last  of  the  inverse-moment-matrix  listings  shown 
in  Figure  6-14.  At  the  bottom  of  Figure  6-14,  the  Lagrange  multipliers 
obtained  from  Eq .  (3-19)  on  page  3-9  of  Reference  4  are  listed. 


Some  scaled-variable  results  from  a  constrained  search  step  are  shown  at 
the  top  of  Figure  6-15.  The  increments  to  the  search  variables,  as  given  by 
the  fiist  of  Eqs.  (3-20)  on  page  3-9  of  Reference  4,  are  listed  first,  fol¬ 
lowed  by  the  contributions  to  these  increments  of  each  of  the  two  terms  on  the 
right  side  of  this  equation. 

The  line  below  these  arrays  is  associated  with  constraint-plane  and 
local-minimum  tests.  The  first  number,  LAM(J) ,  is  the  algebraically  largest 
nonzero  Lagrange  multiplier,  and  the  second,  LMK(J,J)  is  the  corresponding 
diagonal  element  of  the  inverse- moment  matrix.  The  third  number  is  defined  by 
BETA  =  L  *  LAM(J)  *  |l,MK(J,J)|.  The  fourth  number,  PNORM,  is  the  square  root 
of  the  sum  of  the  squares  of  the  search-variable  increments.  A  value  of  PNORM 
less  than  or  equal  to  BETA  indicates  that  the  constraint  plane  corresponding 
to  the  index,  J ,  in  the  above  equation  for  BETA  should  be  dropped.  The 
scaled-variable  increments  are  then  recomputed  with  the  correspondingly 
revised  active- constraint  matrix,  G.  The  values  of  BETA  and  PNORM  also  are 
used  in  the  first  and  second  parts,  respectively,  of  the  three-part  test  shown 
on  page  3-11  of  Reference  4  t  determine  if  a  local  constrained  minimum  has 
been  reached.  (Since  the  factor  h  and  the  inverse-moment  element  in  the 
expression  for  BETA  are  together  of  order  of  magnitude  1,  and  since  the 
Lagrange  muitipiier  used  in  BETA  is  the  largest  of  the  Lagrange  multipliers, 
the  tests  or  the  Lagrange  multipliers  shown  in  Eq.  (3-27)  of  Reference  4  are 
,  pprox  ’  i.ta  ce  1  y  equivalent  to  a  test  ot  BETA). 

The  next,  qu  unity,  PG ,  in  Figure  6-15  is  the  scalar  product  of  the  scaled 
s'  i rch-va rub  1 o  increments  listed  at  the  top  of  the  figure  and  the  scaled 
f lurter-spc  derivatives  in  Figure  6-13.  For  a  positive  or  zero  value  of 

Lis  quantity.  th<>  move  represented  by  these  increments  is  ineffective,  i.e., 

’be  t  lute*,  r  speed  would  be  increased  or  remain  constant.  In  this  situation,  a 
new  seat'a  direction  is  computed  based  on  the  projection  of  the  gradient  on  ^ 

•  he  active  constraint  planes.  . 
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igure  6-14  -  Typical  Listing  of  Search-Module  Output  Pertaining  to 

Introduction  of  Constraints  (results  shown  are  for  different 
analysis  cycle  than  those  in  other  figures  containing  search- 
module  output) . 


scale  factors  provided  as  input  data  via  Items  81  and  8K  of  the  vibration- 
analysis  module.  As  mentioned  previously  in  the  discussion  of  Figure  6-12, 
this  and  other  scaled  arrays  contain  elements  only  for  the  store  parameters 
for  which  the  input  scale  factors  are  nonzero.  The  sum  of  the  squares  of  the 
scaled-derivative  elements  is  shown  at  the  bottom  of  Figure  6-13. 

The  array  immediately  after  the  new  inverse  Hessian  in  Figure  6-13  is  the 
product  of  the  latter  array  and  the  scaled-derivative  array.  The  result  is 
the  direction  array  given  by  Eq .  (3-6)  on  page  3-3  of  Reference  4  (without  the 
minus  sign).  The  last  array  illustrated  in  Figure  6-13  defines  the  search 
point  for  which  the  analyses  and  derivative  computations  have  just  been 
completed.  The  immediately  preceding  array  gives  similar  information  for  the 
previous  analysis  cycle.  E\,i.  the  first  cycle,  the  previous  and  present  arrays 
are  identical. 

The  next  group  of  search-module  output  pertains  to  the  introduction  of 
constraints,  and  is  illustrated  in  Figure  6-14.  For  a  first  pass  through  the 
search  module,  or  if  the  number  of  active  constraints  was  previously  zero,  the 
position  of  the  present  search  point  relative  to  the  constraint  planes  is 
calculated  by  substituting  the  scaled  variables  into  the  expression  on  the 
Left  side  of  Eq.  (3-1)  on  page  3-1  of  Reference  4.  The  result  of  this 
calculation  is  the  first  array  illustrated  in  Figure  6-14.  A  moderate-size 
negative  value  for  an  element  indicates  that  the  present  search  point  is 
inside  the  search  space  with  respect  to  the  constraint  plane  having  an  index 
corresponding  to  that  element.  An  element  value  greater  than  -0.001  is 
roiis f.dere  1  within  the  program  to  denote  that  the  search  point  is  on  or  outside 
the  constrair.t  plane  corresponding  to  that  element.  In  the  sample  output 
shown,  constraints  2,  9,  and  10  are  active. 

/.  projection  of  the  search  direction  along  the  active  constraints  is 
performed  as  discussed  in  Subsection  3.3  of  Reference  4.  As  each  successive 
constraint  is  introduced  (except  the  first),  a  listing  is  given  of  the 
inverse -moment  matrix  (see  page  3-9  of  Reference  4)  corresponding  to  all  the 
constraints  that  have  been  introduced  to  that  stage  of  the  calculations.  If 
no  new  constraints  are  being  introduced  in  a  particular  analysis  cycle,  the 
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The  next  two  arrays  in  Figure  6-12  are  the  changes,  relative  to  the 
previous  search  point,  in  the  search  variables  and  in  the  partial  derivatives 
of  the  flutter  speed  with  respect  to  the  search  variables.  For  these,  as  well 
as  the  subsequent  arrays  in  Figure  6-12,  the  size  of  the  arrays  is  determined 
by  the  number  of  nonzero  scale  factors  entered  in  Item  HI  of  the  search-module 
input  data. 

Next,  the  inverse -Hessian  matrix,  l.e.,  the  inverse  of  the  matrix  of 
second  partial  derivatives  of  the  flutter  speed  with  respect  to  the  search 
variables,  is  listed  for  the  previous  search  point.  Since  two  search  steps 
are  required  to  develop  the  inverse  Hessian  by  the  procedure  used  in  ESP,  the 
default  identity  matrix  will  be  listed  if  the  current  search  point  was  reached 
following  the  first  search  step. 

Below  the  old  inverse -Hess inn  matrix  are  several  intermediate  results 
obtained  during  the  calculation  of  the  updated  version  of  this  matrix.  The 
approach  used  to  calculate  the  inverse  Hessian  is  discussed  in  Subsection  3.2 
of  Reference  -t ,  and  the  final  expressions  used  are  given  as  Eqs.  (3-12)  on 
page  3-6  of  that  reference.  First,  the  values  of  the  second  term  in  the 
second  of  Eqs.  (3-12)  is  listed,  and  then  the  sum  of  both  terms  in  that 
equation  is  given.  The  last  two  quantities  are  the  value  of  c  from  the  third 

■  f  Fqs.  (3-12),  and  the  scalai  product  of  the  previous  Z  array  (y  in  Eqs. 
(3-12))  with  itself.  The  latter  quantity  is  used  in  a  comparison  test  with  c 
to  determine  if  the  inverse-Hessinn  increment  given  by  the  second  term  in  the 
i  fast  of  Eqs,  (3-12)  will  be  unacceptably  large.  It  it  Is,  Z2  Is  used  in 
place  of  i:  in  this  term.  Messages  are  printed  for  this  and  other  tests  if  a 

,<!  ■  i  lent  ion  ot  the  first  of  Eqs.  (3-12)  is  appropiiate. 

The  new  iuve  i.'se-Hess  1  ari  array,  generally  obtained  according  to  the  first 
qs  1  3  -IV),  i '■  shown  at  the  top  of  Figure  6-13.  For  the  first  entry  into 
the  seal  ^h  module ,  this  will  be  the  first  output  after  that  illustrated  In 
igure  6-1 ! ,  and  the  new  inverse  Hessian  will  be  the  identity  matrix.  The 

■  ;  i  pH- Her ’ vat  ive  array,  which  is  the  third  group  of  output  in  Figure  6-13,  is 

■  'f -lined  ’>*,  combining  the  derivatives  illustrated  in  Figure  6-11  with  the 
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0. 

0. 

1. OOOOOE+OO 

0. 

0. 

J  ♦ 

0* 

a 

0. 

0. 

1. OOOOOE+OO 

0. 

J  . 

c. 

0. 

0. 

0. 

1. OOOOOE+OO 

0  * 
o. 

0. 

0. 

0. 

0. 

1  OOOOOE+OO 

H0L0*(DELTA 

1  GRAD) 

3.7492RE-02 

oo0905E-02 

3.311 72E-02 

8 . 50673E-02 

-1 .21997E-01 

3 . 1 7798E-02 

. -  OFLTAX-HOLDG 

s . 1  728  4E -03  - 

7 . 53373E-02 

6 , 99686E-02 

-1.01 752E-02 

1 . 21997E-01 

-8.1 7798E-02 

:k  and  Z2  = 

- . 33017E 

-01  . 32322E-01 

i\\iro  A— 1 2  -  Typical  Listing  of  Search-Module  Output  Defining  Constraint 

Status,  Search  Progress,  and  Intermediate  Results  Required  for 
Updating  the  Inverse-Hessian  Matrix  (Obtained  Only  After 
Second  or  Subsequent  Analysis  Cycle). 
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degrees  of  freedom  (see  Eq .  (4-1),  page  4-1,  Reference  4).  The  vector  com¬ 
ponents  are  listed  as  pairs  of  real  and  imaginary  parts  in  an  order  corres¬ 
ponding  to  the  store  degree-of-f reedom  numbers  (cf.  Figure  6-4  and  associated 
discussion).  The  units  of  the  rotation  degrees  of  freedom  are  rad. /ft., 
assuming  a  normalization  with  respect  to  displacement.  The  quantity  CSCt  is 
defined  by  Cq .  (9./),  page  107,  ,-f  Reference  1,  Volume  1;  it  is  a  repeat  of 

the  last  listed  item  be  foie  the  print-plots  in  the  flutter-analysis  module. 

The  next  group  of  output  n  he  described  (see  Figure  6-12)  appeals  only 
alter  a  search  step  has  been  taken ,  i.e.  ,  in  a  second  or  subsequent  entry  to 
the  search  model:;.  At  the  top  of  the  figure,  the  status  of  the  search  with 
respect  to  the  various  constraint  planes  is  given.  The  numbering  system 
corresponds  to  the  order  in  which  constraint  planes  are  defined  by  means  of 
either  the  groups  of  constraint-equation  parameters  entered  in  input-data 
items  oP  and  8b  j T  T0C1  -  1  ,  or  the  points  on  the  constraint  planes  entered 
,u  lift,!  8.  (to:  i  1 01.  -  It,  It-!:  8.0?/  { for  I  TOC  *  2),  and/or  Item  8GG  (for  HOC 
---  J)  , f  l  IOC  t  <J. 

The  next  : ten:  in  figure  6-12  is  the  inverse-moment  matrix  defined  near 
the  top  of  page  T-c>  of  Reference  4.  The  version  of  this  matrix  that  is  listed 
■a  this  figure  was  computed  in  che  previous  search  step,  and  the  values  in 
figure  o-12  are  a  repetition  of  the  values  provided  toward  the  end  of  the 
oct put  for  that  step  (see  Figure  6-i  4  and  associated  discussion  below).  The 
size  of  die  inverse-moment  matrix,  3x3  for  the  case  shown  in  Figure  6-12, 
corresponds  to  the  number  of  active  constraints.  This  matrix,  as  well  as  all 
•>ubp '••queri'"  .-•'•arch  -medvie  printed  output  except  for  flutter  speeds  and  the  new 
a.:  parameters,  are  given  in  terms  of  scaled  variables. 

fee  I  1  utter -speed  change  1  is  ted  next  is  with  respect  to  the  value  for  the 
.!  f.  ions  ana lysis  cycle.  This  quantity  is  listed  for  the  first  entry  to  the 
■■I  ■  -  .a.  T:  if/  /  .3  vo  1 1  as  for  the  second  and  subsequent  entries,  but  for  the 
iii.,:  t-ntrv ,  the  p  >  .•need  ..suit  is  m  :n  i  ugless  in  that  it  is  with  respect  to 
:  t  f  1  s  1  *  v  ;  i  i  1  nigh  ini  tin!  units'  specified  ’or  VNKW  in  1 1  eni  RK  ot  the 
i  brat  :  o-.ina  1  vs  ’  s-.t.  .'du  1  ■?.  input  data. 


tv  -  i  6 


The  reader  is  reminded  that,  as  indicated  in  Items  8A,  8J,  and  8K  in 
Subsection  5.2,2  of  this  report,  a  restart  run  requires  certain  changes  to  the 
primary  input-data  file  in  addition  to  the  use  of  the  restart  file  as  an  input 


via  TAPE47. 
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APPENDIX  A 


NASTRAN/ESP  INTERFACE 


A. I  -  COSMIC  NASTRAN 

As  was  discussed  briefly  in  Subsection  5.1,  one  form  in  which  dynamics- 
model  matrices  may  be  read  in  ESP  is  that  provided  by  the  0UTPUT2  routine  in 
the  COSMIC  version  of  NASTRAN  (see  Reference  8,  pages  5.5-24  through  5.5-27). 
In  this  interface  as  presently  implemented,  any  matrices  that  are  to  be 
obtained  from  COSMIC  NASTRAN  must  be  written  via  the  same  output  unit.  In  the 
control-statement  sequence,  the  local-file  designator  for  this  unit  is  UT1, 
and,  in  the  required  additional  DMAP  statements,  the  unit  designator  should  be 
set.  to  11. 

It  is  noted  parenthetically  that  an  interface  in  which  matrices  are 
written  to  separate  files  (via  different  output  units)  would  have  provided  the 
user  with  a  greater  freedom  of  choice  in  selecting  combinations  of  ESP  input 
matrices.  However,  at  least  for  the  CDC  April  1984  COSMIC  NASTRAN  release, 
for  which  the  present  Interface  was  developed,  the  code  as  delivered  provides 
for  only  one  output  unit  in  conjunction  with  0UTPUT2.  A  Fortran  modification 
could  be  introduced  to  obtain  the  more  general  capability  described  in 
|  Reference  8.  However,  it  was  judged  that  most  installations  would  probably  be 

■isfng  the  program  as  received  from  COSMIC,  and  that  users  at  these 
:  a- r,-i.  1  lations  would  prefer  to  restrict  their  0UTPUT2  usage  to  a  single  unit 
ither  than  change,  or  request  a  change  to,  their  original  NASTRAN  source 


'  is  noted  also  that  difficulty  was  encountered  in  obtaining  results 
.  w;'Tpi'T?  using  release  17.7  of  COSMIC  NASTRAN.  Thus  far,  the  procedure 
described  herein  has  been  used  successfully  only  with  the  April  1984  release. 

An  illustrative  control -statement  sequence  for  executing  the  April  1984 
.■el  ease  rl  COSMIC  NASTRAN  on  the  NADC  Central  Computer  System  is  shown  in 
; igore  A- L .  In  addition  to  containing  the  time  parameter,  the  JOB  card 


NASTRAN»TXXXX»STNOS. 
ACCOUNT  »  UUUUUU  *  PPPPPP . 
GET » DDDDDDD . 

DEFINE, UT1=KNDM. 

GET  » NAST484/UN=SYSTEM . 
BEGIN, »NA3T484, 

RFL » 160000 « 

REDUCE ( - ) . 

L INK1 »  DDDDDDD . 


Typical  Control-Card  Sequence  for  COSMIC  NASTRAN  Run  t 
Prepare  0UTPUT2  File  for  ESP. 


specifies  execution  via  the  NOS  operating  system.  Execution  is  initiated  via 
the  LLNK1  statement,  in  which  the  second  term  is  the  input  data  file. 

The  permanent  file  name  KMDM,  shown  in  conjunction  with  local  file  UT1  in 
rigure  A-l ,  is  intended  as  a  reminder  of  the  required  matrix  order  in  the 
t!:TPl'T2  file.  As  will  be  seen  from  the  subsequent  discussion  of  the  addi- 
r tonal  DMAT  statements,  each  letter  in  the  OUTPUT2  permanent-file  name  is  the 
!  List  letter  of  a  DMAP  name  that  has  been  assigned  herein  to  one  of  the  four 
i  it  rices  used  by  ESP: 

o  K  -  flexibility  matrix 

o  M  -  dynamic  mass  matrix 

.<  D  -  rigid-body-displacement  matrix 

o  M  -  plug  mass  matrix 

io'ii'  however  that,  as  indicated  in  Subsection  5.1,  not  all  ESP  input  matrices 
v:-  t  be  read  from  the  same  source,  and  therefore  the  0UTPUT2  file  should  not 
.M'Cessarily  include  all  four  of  the  matrices  listed. 

Figure  A-2  shows  a  representative  set  of  DMAP  modifications  to  Rigid 
Format  3  (Normal  Mode  Analysis)  to  obtain  a  file  of  the  desired  ESP  matrices 
via  0UTPIJT2 .  The  statements  to  be  added  to  the  Executive  Control  Deck  are 
.■how  in  part  (a)  of  the  figure.  Part  (b)  shows  these  statements  in  the 
•  or. text  of  the  original  DMAP  sequence,  so  that  the  relationship  of  the  alters 
'  the  rigid  format  in  which  they  are  inserted  can  be  seen  without  having  to 
Tvtvr  either  to  the  printed  DMAP  sequence  in  Reference  8  (pages  3.4-1  through 
or  to  a  listing  obtained  from  a  NASTRAN  execution  using  a  DIAG  14 
..srcutive  Control  card.  For  a  more  complete  understanding  of  the  relationship 
:  t'.vcc  the  alters  and  the  original  DMAP  sequence,  the  reader  may  wish  to 
v i t  the  description  of  the  DMAP  operations  in  Rigid  Format  3  given  on 
,  <  4-7  through  3.4-1]  of  Reference  8. 

A  key  consideration  in  determining  the  DMAP  alters  for  providing  the  ESP 
:,  i{  rices  is  the  set  of  units  used  in  the  NASTRAN  analysis.  Independent  of  the 
i'-urce  of  the  needed  matrices,  the  units  must  be  as  specified  in  Subsection 
l  For  the  case  illustrated  in  Figure  A-2,  the  Bulk  Data  was  in  the  same 
,  '  as  those  used  in  ESP,  and  a  PARAM  WTMASS  card  had  been  included  to 
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ALTER  66 

SOLVE 

MATPRN 

PARAMR 

PARAMR 

ADD 

ADD 

MATPRN 
ALTER  68 
MATPRN 
OUTPUT2 
0UTPUT2 


KLL»/KLLI/  ♦ 

NLLI t *, >//  % 

//*DIV*/V»N»RECUTM/1 .0/C» Y»UTMASS  % 
//♦COMPLEX#/ /V  »N»RECUTM/0.0/V>N»RECWTMC  * 
MLL»/MLLW/VfN*RECUTMC  $ 

MRR f /MRRU/V >N>RECWTMC  % 

MLLU » MRRW  *  »  » //  « 

DM »  »  »  » //  $ 

KLLI»MLLW>DM»MRRU»//-1/11  % 
t,,?//- 9/11  « 


(a)  Additions  to  Executive  Control  Deck 


62 

SMP1 

USET»KFF» , »/G0»KAA,K00»L00»  »  >  »  »  % 

63 

SMP2 

USET»GO»MFF/MAA  $ 

64 

LABEL  • 

LBL5  $ 

65 

COND 

LBL6  t REACT  * 

66 

RBM61 

USET»  KAA.MAA/KLL  »KLR»KRR»  MLL  »MLR»MRR  % 

66 

SOLVE 

KLL » /KLLI/  $ 

66 

MATPRN 

KLL I »»»»//  $ 

66 

PARAMR 

//*DIV*/V,N»RECUTM/1 , 0/C r Y , UTMASS  $ 

66 

PARAMR 

//♦COMPLEX*//" »N»RECWTM/0.0/V»N»RECWTMC  * 

66 

ADD 

MLL»/MLLU/V>N,RECUTMC  * 

66 

ADD 

MRR »/MRRW/V»N»RECWTMC  * 

66 

MATPRN 

MLLU  »  MRRW »  >  » //  $ 

67 

RBMG2 

KLL/LLL  * 

68 

RBM63 

LLL»KLR»KRR/DM  % 

68 

MATPRN 

DM »»»»//  * 

68 

0UTPUT2 

KLLI>MLLW»DM»MRRW>//-1/11  * 

68 

OUTPUT2 

t  t  t  t //-9/1 1  $ 

69 

RBMG4 

DM»MLL  » MLR r MRR /MR  * 

70 

LABEL 

LBL6  $ 

(b)  Modified  Rigid-Format-3  Listing 


Figure  A-2  -  Illustrative  Modifications  to  COSMIC  NASTRAN  Rigid  Format  3 
(Normal  Mode  Analysis)  to  Obtain  Desired  ESP  Matrices 
via  0UTPUT2 


convert  the  weight  data  to  mass  units  within  NASTRAN.  To  return  to  the  weight 
units,  the  reciprocal  of  the  WTMASS  parameter  (RECWTM)  was  computed,  and  this, 
when  converted  to  the  required  complex  form,  was  used  to  multiply  the  dynamic 
mass  matrix  (MIL)  and  the  plug  mass  matrix  (MRR) . 

To  be  consistent  with  the  procedure  for  reading  the  COSMIC  NASTRAN 
matrices  in  ESP,  the  first  parameter  in  the  0UTPUT2  calling  sequence  (see 
Reference  8,  page  5.5-24)  must  be  -1  when  writing  the  desired  matrices.  The 
second  parameter,  which  is  the  output  unit  number,  will,  as  noted  previously, 
usually  be  11  for  CDC  COSMIC  NASTRAN.  The  third  parameter  is  not  needed,  and 
is  therefore  omitted  in  the  example  shown  in  the  figure.  Following  the  use  of 
01TPUT2  to  write  the  matrices,  a  second  0UTPUT2  statement  is  included  (accord¬ 
ing  to  instructions  in  Reference  8)  to  write  an  end-of-file  mark. 

As  a  final  comment  on  obtaining  ESP  matrices  from  COSMIC  NASTRAN,  and  MSC 
NAS  i’RAN  as  well,  it  is  noted  that  the  dynamics  idealization  used  must  be 
chosen  so  as  to  be  consistent  with  the  ESP  interpolation  procedure.  Specifi¬ 
cally,  the  dynamics  grid  points  must  lie  on  a  set  of  spanwise-oriented  lines 
as  illustrated  in  Reference  1,  Volume  I,  page  86. 

To  complete  the  discussion  of  the  interface  between  COSMIC  NASTRAN  and 
ISP,  a  sample  of  the  printed  output  that  is  obtained  from  ESP  when  matrices 
are  read  from  an  0UTPUT2  file  is  shown  in  Figure  A-3.  The  printing  of  the 
flexibility  and  dynamic  mass  matrices  is  controlled  by  KLUEV(7)  and  KLUEV(8), 
respectively  (see  Subsection  5.2.2);  the  rigid-body-displacement  matrix  and 
the  plug  mass  matrix  are  always  printed.  The  matrix  name  shown  in  the  ESP 
listing  is  taken  from  the  name  used  in  the  0UTPUT2  DMAP  statement. 

A. 2  -  MFC  NASTRAN 

The  interface  between  the  MacNeal-Schwendler  Corporation  (MSC)  version  of 
NASTRAN  and  ESP  parallels  that  developed  for  COSMIC  NASTRAN  except  for  a  few 
variations.  First,  to  provide  greater  freedom  of  choice  in  selecting  combina¬ 
tions  of  FSP  input  matrices,  the  four  desired  matrices  are  assumed  to  be 
written  by  MSC  NASTRAN  via  0UTPUT4  to  three  different  output  units  rather  than 


r  j  -t.  ro  <_n 


READING  COSMIC  NASTRAN  MATRIX  DM 
NUMBER  OF  COLUMNS  =  3 

NUMBER  OF  ROUS  =  130 


READING  COLUMN  1 
NUMBER  OF  ROWS  READ  =  130 

1 • OOOOE+OO  3 • 544 IE- 1 4  -5.4770E-13  l.OOOOE+OO  3.4475E-14 
, 2374E- 1 5  -2 . 3313E-1 3  l.OOOOE+OO  9.4279E-16  -4.5516E-14 
. 2774E- 10  l.OOOOE+OO  -1.2897E-11  2.2778E-10  l.OOOOE+OO 

.OOOOE+OO  - 1 . 2897E-1 1  2.2781E-10  4.1930E-07  -5.6440E-09 

. 64 1 1 E-l 2  2 . 9308E-07  -2.2524E-09  -3.0687E-10  -9.7966E-14 

• 6587E- 1 0  -3 • 2529E-1 1  4.8359E-10  -2.1033E-10  -6.1954E-11 

- 1 . 4661E-1 1  5.9053E-10  -4.5631E-11  -4.9853E-12  6.0109E-10  ••• 

-1 . 2362E-10  7 • 1660E-1 1  1.6474E-09  -5.2270E-02  -1.5679E-12 

-3  •  8408E- 1 1  2 . 6592E-09  -4.2439E-10  -3.8408E-U  1.1875E-08 

-4 . 6407E-02  -5.5957E-10  -4.6407E-02  5.2812E-10  -4.6407E-02 

-4 . 6407E-02  -1.5956E-10  -3.30S7E-02  1.2539E-09  -3.3057E-02 

8 . 2763E-09  -5.4285E-11  6.3358E-09  -5.7732E-11  4.5312E-09 

3 . 1 170E-09  -1 . 2079E- 10  3.4353E-09  -1.7612E-10  4.6702E-09 


READING  COLUMN  2 
NUMBER  OF  ROWS  READ  =  130 

-6 • 3810E+00  l.OOOOE+OO  1.5607E-12  -3.7010E+00  l.OOOOE+OO 

l.OOOOE+OO  5 . 1404E-13  1.6552E-12  l.OOOOE+OO  9.2646E-14 

1 . 6479E-09  -1.5791E-07  l.OOOOE+OO  1.6479E-09  -2.0545E-07 

-2.8628E-07  l.OOOOE+OO  1.6479E-09  3.1751E-06  -2.1575E+02 

l.OOOOE  +  OO  2  •  1 656E-06  -1.8100E  +  02  8.8518E-08  l.OOOOE  +  OO  ••• 

1 . 1644E-07  l.OOOOE+OO  -1.2600E+02  1.0246E-07  l.OOOOE+OO 

l.OOOOE+OO  -6 . 6430E+01  4.6917E-08  l.OOOOE+OO  -4.9675E+01 

1 . 9330E-08  9 . 9863E-0 1  -5.2340E-02  -8.8280E+01  -3.0605E-08 

l.OOOOE+OO  -2 . 2900E+02  -9.8517E-08  l.OOOOE+OO  -2.2900E+02 


Figure  A-3  -  Typical  ESP  Listing  Obtained  When  Reading  COSMIC  NASTRAN 
0UTPUT2  File 


t  —  w—i  —  v-'v.'-'i  ii p  t  my 


only  one  as  with  COSMIC  NASTRAN.  Three  units  should  be  used,  rather  than 
four,  because  the  plug  mass  matrix  is  considered  to  be  appended  to  the  file 
containing  the  dynamic  mass  matrix. 

A  second  difference  is  that  the  ESP  interface  with  MSC  NASTRAN  has  been 
written  to  take  advantage  of  the  two  matrix  sparseness  options  that  are 
available  in  the  0UTPUT4  routine  (see  Reference  9,  pages  5.4-91  through 
5.4-92a).  The  usually  sparse  dynamic  mass  matrix  is  assumed  to  be  transferred 
using  the  sparse  option,  and  the  transfer  of  the  other  three  matrices  is  done 
via  the  nonsparse  option. 

A  typical  set  of  MSC-NASTRAN  DMAP  alter  statements  for  Rigid  Format  3  in 
CUC  Release  63,  April  1983,  is  shown  in  Figure  A-4,  part  (a).  To  show  the 
context  in  which  these  statements  are  used,  portions  of  the  resulting  modified 
DMAP  listing  are  given  in  part  (b)  of  this  figure.  A  comparison  of  the 
0HTPUT4  statements  contained  in  the  figure  with  the  0UTPUT4  description  given 
in  Reference  9  illustrates  the  points  made  in  the  previous  paragraph. 

Figure  A-5  shows  the  form  of  the  printed  output  from  ESP  when  reading  a 
matrix  from  an  0UTPUT4  file.  The  "type"  of  the  MSC  NASTRAN  matrices  (see 
Reference  9,  page  5.4-2a)  is  required  by  ESP  to  be  real,  single-precision, 
i n d ,  to  provide  a  check  that  this  rule  has  been  followed,  the  numerical  index 
for  the  type  is  listed  from  the  first  record  of  each  0UTPUT4  file.  The  header 
information  for  each  column,  which  is  obtained  from  the  beginning  of  the 
OUTPUT 4  record  for  that  column,  provides  additional  information  for  checking 
purposes  related  to  the  data-condensation  features  in  0UTPUT4. 
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ALTER  181 

SOLVE  KLL»/KLLI/  * 

MATPRN  KLLI//  $ 

0UTPUT4  KLLI//-3/11  * 

ALTER  212 
MATPRN  DM//  * 

0UTPUT4  DM//-3/12  * 

ALTER  250 

PARAMR  //DIV/V*N»RECWTM/1 . O/WTMASS  $ 

PARAMR  //COMPLEX/ /RE CWTM/O  »0/V»N»RECWTMC  $ 

ADD  MLL » /MLLW/V  »N»RECWTMC  t 

ADD  MRR»/MRRW/V»N»RECWTMC  * 

MATPRN  MLLU »  MRRW//  t 
0UTPUT4  MLLW//-1/-13  $ 

OUTPUT  4  MRRW//-2/13  $ 


(a)  Additions  to  Executive  Control  Deck 
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UPARTN 
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• 

• 

USET  »KTT/KLLt  » KLR » KRR/T/L/R  $ 

181 

SOLVE 

KLL » /KLLI/  * 
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MATPRN 

KLLI//  i 
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0UTPUT4 

KLLI//-3/11  * 

212 

RBMG3 
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LLL*  »KLR»KRR/DM  * 

212 

MATPRN 

DM//  * 
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OUTPUT  4 

DM//-3/1 2  * 
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UPARTN 
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USET  »MTT/MLL»  » MLR » MRR/T/L/R 

$ 
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PARAMR 

//DIV/V>N»RECWTM/1 .O/WTMASS 
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PARAMR 

//COMPLEX //RE CWTM/O .  0/V  » N  f  RECWTMC  * 
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ADD 
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250 
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0UTPUT4 

MLLW//-1/-13  * 
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OUTPUT  4 

MRRW//-2/1 3  $ 

251 

RBMG4 

DM,MLL»MLR» MRR/MR  $ 

• 

• 

(b) 

• 

Modified  Rigid-Format-3  Listing 

Figure  A-4  -  Illustrative  Modifications  to  MSC  NASTRAN  Rigid  Format  3 
(Normal  Mode  Analysis)  to  Obtain  Desired  ESP  Matrices 
via  OUTPUT 4 
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Figure  A-5  -  Typical  ESP  Listing  Obtained  When  Reading  MSC  NASTRAN  0UTPUT4  File 


APPENDIX  B 


IMPLEMENTATION  OF  CAPABILITY  FOR  RIGID-BODY  MODES 

Since  a  capability  for  representing  free-free  modes  was  contained  in  the 
FASTOP  code  from  which  the  original  version  of  ESP  was  developed,  a  major 
portion  of  the  capability  for  including  rigid-body  modes  in  the  output  from 
the  vibration-analysis  module  could  be  implemented  by  simply  obtaining  the 
mode  shapes  themselves  from  the  transformation  matrix  from  relative  to 
absolute  coordinates. 

With  the  mode  shapes  being  available,  the  rigid-body  generalized-mass 
terms  could  be  readily  calculated  as  an  integral  part  of  the  calculation  for 
the  the  flexible-mode  generalized  masses.  The  rigid-body  portion  of  the 
generalized-mass  matrix  would  generally  contain  off-diagonal  terms,  however, 
since  the  ignorable  coordinates,  or  "plug"  degrees  of  freedom  (see  Reference 
1,  Volume  I,  Section  7,  pages  48  and  49),  would  generally  not  be  coincident 
with  the  airplane  center  of  gravity. 

With  the  introduction  of  a  capability  for  user  specification  of  zero- 
::  1 r speed  rigid-body  frequencies  (see  page  3-2),  a  procedure  for  calculating 
rigid-body  elements  in  the  generalized-stiffness  matrix  was  required.  It  was 
; nd ge J  most  desirable  for  the  user-specified  rigid-body  frequencies  to  apply 
to  uncoupled  inodes,  i.e.,  to  modes  in  which  the  rotations  are  about  the 
!: plane  tenter  of  gravity.  With  this  approach,  the  effective  rigid-body 
-t-r :  would  be  located  at  the  center  of  gravity,  and  the  general  i  zed-st  1  f  f- 

o  -  si.  matrix,  in  terms  of  modes  defined  by  the  plug  degrees  of  freedom,  would 
■  i  i  '.i'.r  r  ~  the  generalized-mass  matrix  in  that  it  too  would  contain  off- 
f  <■  rrr.p  . 


In  tie  discussion  that  follows,  the  elements  of  the  general ized- 
.  if ‘ness  matrix  are  developed  in  terms  of  the  elements  of  the  generalized- 
r-'s-;  matrix  and  the  specified  natural  frequencies  for  a  rigid-body  system  with 
.x  translation  mode  and  one  rotation  mode.  Thereafter,  expressions  are  given 
far  the  more  general  modal  combinations  associated  with  either  symmetric 
.  i  hi,:  i  t  ud  inn] )  or  antisymmetric  (lateral)  motions. 


If  Che  system  shown  above  has  mass  m  and  moment  of  inertia  I  about  the 
center  of  gravity,  the  equations  of  motion  may  be  written  as 


m  md 

h  ] 

I 

T3 

rC 

\ _ 

h 

)  + 

2 

2 

m  I  +  md “ 

9  ) 

A4  +  khd  . 

0 

0  .  (B-l) 


From  inspection  of  the  figure,  or  from  a  solutio .  of  the  equations  of  motion, 
the  natural  frequencies  of  this  system  are 


,  ,1/2  ,,  /T,  1/2 
^h  =  (^/m)  .  =  (kg/*) 


(  B—  2  ) 


To  develop  the  desired  relationships  between  the  generalized-stiffness 
matrix  elements  and  the  generalized-mass  matrix  elements  and  the  natural 
frequencies,  F  q .  (3-1)  is  first  rewritten  as  follows: 


=  0  .  (3-3) 


By  comparing  Eqs.  (B-l)  and  B-3),  it  it  seen  that  d  may  he  written  as 


AD-A152  268  ESP  (EXTERNAL-STORES  PROGRAM)  -  A  PILOT  COMPUTER  2/2 

PROGRAM  FOR  DETERMINING.  OJ)  GRUMMAN  AEROSPACE  CORP 
BETHPAGE  NV  J  B  SMEDFJELD  FEB  85  ADCR-85-l-VOL-i 
UNCLASSIFIED  N000i9-8i-C-0395  F/G  9/2  NL 
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Now,  by  using  all  of  Che  above  equations,  three  of  the  generalized- 
stif fness/generalized-mass/natural-frequency  relationships  may  be  written 


immediately  as 


shh  '  "hh  “h  • 
'  "he  %  ■ 


M0h  ^h  *  ^10  * 


The  expression  for  '<  requires  some  algebraic  manipulation.  First,  an 

U  0 

expression  for  I  is  formed  by  comparing  Eqs.  (B-l)  and  (B-3)  and  substituting 
Eq.  ( B-4)  : 


I  =  M„ 


Substituting  this  result  into  the  second  of  Eqs.  (B-2)  yields 


■  ("••  -  i)< 


Finally,  after  equating  stiffness  elements  in  Eqs.  (B-l)  and  (B-3),  and 
::ii icing  Eqs.  (B-4),  (B-6) ,  and  (B-9) ,  the  desired  expression  for  K00  is 


ir  jined  : 


K0o  ~ 


M00  “0  + 


is  \  2  /  ^0  \ 

J  "e  +  VW 

(i)  ('  -  4 


(B-10 


For  the  general  symmetric  case  consisting  of  fore-and-aft  translation, 
vertical  translation,  and  pitch,  which  will  be  denoted  here  by  the  subscripts 
i,  2,  and  3,  respectively ,  the  generalized  stiffness  elements  are 


For  the  general  antisymmetric  case  consisting  of  lateral  translation, 
roll,  and  yaw,  denoted  by  the  subscripts  1,  2,  and  3,  respectively,  the 
generalized  stiffness  elements  are: 

Kn  ‘  Mn  “?  •  (1 


mI2  ^  . 


M13  “l  ' 


M23  “1  • 


K22  “  M22  u2  + 


M33  w3  + 


(£)(- :  -  -0  • 
©(-■-4 


The  calculation  of  the  rigid-body  generalized  stiffness  elements  in  ESP 
is  performed  via  Eqs.  (B-ll)  through  (B— 22) . 


APPENDIX  C 


i 

DEFINITION  OF  CONSTRAINT  PLANES 

For  the  results  obtained  from  ESP  to  have  the  most  practical  sipnif icance, 

I  It  Is  desirable  that  a  minimum-flutter-speed  point  in  the  store-parameter  search 

space  be  close  to  an  actual  store  configuration.  This  in  turn  makes  it 
desirable  to  define  the  search  space  in  such  a  manner  that  regions  with  no  or 
few  store  points  are  minimized.  One  step  that  can  be  taken  toward  meeting  this 
objective  is  to  allow  isolated  extreme  store-parameter  combinations  to  fall 
outside  the  search  space,  and  then  analyze  these  points  individually.  Beyond 
that,  however,  the  degree  of  refinement  of  the  search-space  boundaries  will 
usually  be  dependent  primarily  on  the  amount  of  effort  devoted  to  the 
const  rain t -plane  definition. 

O.i  -  FULLY  INDEPENDENT  SEARCH  VARIABLES 

The  simplest  approach  to  the  definition  of  constraint  planes  is  to 
!  assume  that  the  range  of  the  search  for  each  store  parameter  is  independent 

of  the  value  of  every  other  parameter.  With  this  approach,  only  an  upper  and 
a  lower  bound  for  each  search  variable  need  be  specified,  and  the  search 
space,  when  scaled  according  to  the  suggestions  in  Item  81,  Subsection  5.2.2, 
i  will  then  be  a  unit  multi-dimensional  cube  oriented  such  that  each 

search-variable  axis  is  normal  to  one  pair  of  parallel  cube  faces.  Items  8P , 

8R ,  and  8S,  which  are  the  data  items  associated  with  the  general  method  of 
constraint-plane  definition  in  ESP,  will  then  take  on  the  following  values: 

o  The  number  of  constraint  planes,  MSTAR,  which  determines  the  number 
of  columns  in  the  array  G(J,I),  will  be  twice  the  number  of  search 
variables,  NPARM. 

o  Each  unit  normal  vector,  constituting  one  column  of  length  NPARM  in 
the  array  G(J,I),  will  have  only  one  search-space  component,  J,  that 
is  nonzero.  For  the  upper-bound  constraints,  that  nonzero  component 
will  be  +1.C,  and,  for  the  lower-bound  constraints,  it  will  be  -1.0. 

A  sample  G(J,I)  array  for  fully  independent  search  variables  is  shown 
i ti  F igure  C- 1  (a)  . 
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MSTAR 

(2*NPARM) 


NPARM 

=5 


1.0 


-1.0 


1.0 


-1.0 


1.0 


-1.0 


1.0 


-1.0 


1.0 


-1.0 


(a)  Fully  Independent  Search  Variables  (5  variables  shown) 


MSTAR 

=  (2*4)  +  (1*2) 

G(l,4) 

0(2,4) 

G (3 , 5)  ...  G ( 3 , 8) 

G ( 4 , 5 )  ...  G ( 4 , 8 ) 

1.0  -1.0 


(b)  Partial  Two-Dimensional  Search-Range  Dependency  (1  fully  independent 
variable  and  2  store  stations  with  two-dimensional  dependency 
represented  by  4-sided  polygons) 


NPARM  G ( 1 , 1 )  G( 1 , 2) 

*(2*2)  0(2,1)  G(2 , 3) 

+  1 


NPARM 
*  3 


MSTAR 


G  (  1  , 1 ) 


G( 1 , 3) 

0(2,2)  G(2 , 4) 

G (3 , 3) 


G ( 1 , 5) 

0(3,3) 


G ( 1 ,6) 
0(2,6) 


0(1,7) 
G  ( 2 , 7) 


0(2,8) 

0(3,8) 


(c)  Three-Dimensional  Search-Range  Dependency  (1  store  station  with 
dependency  represented  by  polyhedron  with  8  faces) 


Figure  C - 1  -  Illustrative  Unit-Normal-Vector  Arrays,  0(1,1) 


o  In  Item  8S,  the  scaled  distance,  B,  to  each  constraint,  I,  will  be 
measured  along  the  search-space  axis  corresponding  to  the  single 
nonzero  component  of  the  associated  unit  normal  vector. 

0.2  -  VARIABLES  WITH  TWO-DIMENSIONAL  SEARCH-RANGE  DEPENDENCY 

In  general,  an  effort  toward  some  refinement  of  the  search  space  will  be 
warranted,  and,  if  the  refinement  is  limited  to  two  variables  (e.g.,  mass  and 
pitch  inertia)  at  each  station,  it  is  easily  introduced  in  many  situations.  A 
key  consideration  is  whether  the  range  of  store  parameters  at  one  station  can  be 
assumed  to  be  unaffected  by  the  stores  at  other  stations.  If  so,  the 
constraint  equations  for  each  station  can  be  developed  independently,  and  the 
two  variables  at  each  station  for  which  the  search  ranges  are  to  be  mutually 
dependent  will  constitute  a  two-dimensional  space  such  as  the  schematic  shown 
in  Figure  3-1,  page  3-2,  Reference  4.  Any  other  variables  at  each  station, 
for  which  no  search-space  refinement  is  to  be  introduced,  will  be  independent 
of  each  other  and  also  independent  of  the  two  variables  involved  in  the 
ref inement . 

The  determination  of  the  contributions  to  G(J,I)  and  B(I)  from  each 
two-dimensional  space  may  be  accomplished  using  either  of  two  approaches.  A 
primarily  graphical  approach  involves:  (I)  drawing  an  outward  normal  from  the 
scaled  search  space  at  each  polygon  side,  and  then  determining  its  direction 
cosines  (components  of  the  unit  normal  vector);  and  (2)  measuring  the  scaled 
normal  distance  from  the  origin  to  each  polygon  side  or  extension  thereof. 
Alternatively,  the  scaled  polygon  vertices  may  be  used  to  develop  an  equation 
n-  each  polygon  side;  then  these  equations  may  be  manipulated  into  the  form 
in  Eq .  (3-1),  page  3-1,  Reference  4,  from  which  the  desired  results  are 
i mred i  Italy  available . 

Th"  data  in  Items  8P,  8R,  and  8S  may  now  be  obtained  as  a  combination  of 
the  results  for  the  individual  store  stations: 

o  MSTAR,  which  corresponds  to  the  number  of  columns  in  the  array 
G ( J , I ) ,  will  be  the  sum  of  two  contributions:  (1)  the  sum  of  the 
number  of  sides  associated  with  the  polygons  for  the  two  store 
parameters  at  each  station  with  mutually  dependent  search  ranges;  and 


(2)  twice  the  number  of  additional  mutually  independent  variables  (as 
in  Subsection  C.L). 

o  The  columns  of  G(J,I),  which  are  composite  unit  normal  vectors 
considering  all  search  variables  at  all  store  stations,  will  in 
general  have  two  nonzero  elements  for  those  vectors  associated  with  a 
polygon  side.  Vectors  associated  with  constraints  for  mutually 
independent  variables  will  have  one  nonzero  component  (as  in 
Subsection  C.l).  A  G(J,I)  array  illustrating  this  combination  is 
shown  in  Figure  C- 1(b). 

o  The  array  B ( I )  wilL  be  comprised  in  part  of  elements  equal  to  the 
distances  to  the  polygon  sides  and  in  part  of  elements  similar  to 
those  in  Subsection  C.l. 

C . 3  -  VARIABLES  WITH  THREE-DIMENSIONAL  SEARCH-RANGE  DEPENDENCY 

If  the  refinement  of  a  store-parameter  space  is  to  be  carried  a  step 
further,  such  that  the  search  ranges  of  three  variables  are  considered 
mutually  dependent,  the  previous  discussion  might  lead  to  the  conclusion  that 
three-dimensional  modeling  is  required.  However,  as  is  discussed  and 
illustrated  in  Subsection  7.2  of  Reference  4,  an  alternative  procedure  is  to 
define  the  desired  enclosing  polyhedron  for  a  three-dimensional  space  (e.g., 
mass,  M,  pitch  inertia,  I,  and  longitudinal  center  of  gravity,  x)  by 
constructing  constraint  lines  in  three  two-dimensional  spaces  ((M,I),  (M,r), 
and  ( x , I ) )  successively. 

An  important  advantage  of  this  approach,  other  than  facilitating  the 
definition  of  the  three-dimensional  search  boundary,  is  the  fact  that  each  face 
of  the  resulting  polyhedron  appears  as  a  line  in  at  least  one  of  the  throe 
wo-d  i  mens  i  ona  1  spaces.  Thus,  the  unit  normal  vectors  to  he  provided  in  the 
o(.l, I)  array  will  each  consist  of  no  more  than  two  nonzero  components  which  can 
he  determined  from  two-dimensional  plots  exactly  as  described  in  Subsection  C.2 

The  data  in  Items  8P,  8K,  and  8S  of  Subsection  5.2  will  again  consist  ol 
a  combi  pat  ion  of  the  data  for  individual  stations.  As  in  the  previous 
subsection,  fully  Independent  search  variables  may  be  added  to  those  having 
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dependency,  and,  if  desired,  two-dimensional  dependency  may  be  assumed  at  one 
store  station  and  three-dimensional  dependency  at  another.  For  this  more 
general  situation: 

o  The  number  of  constraint  equations,  MSTAR,  will  be  the  sum  of:  (1) 
the  total  nu.ber  of  polyhedron  faces;  (2)  the  total  number  of 
polygon  sides;  and  (3)  twice  the  total  number  of  mutually 
independent  search  variables. 

o  The  unit  normal  vectors,  which  comprise  the  columns  of  G(J,I),  will 
again  have  at  most  two  nonzero  components.  However,  the  two 
nonzero  components  for  a  polyhedron  face  may  correspond  to  any  two  of 
the  three  search  parameters  associated  with  that  polyhedron,  and, 
therefore,  as  illustrated  in  Figure  C-l(c),  the  two  nonzero  elements 
in  a  column  of  G(J,I)  will  not  necessarily  be  contiguous. 

o  Since  the  polyhedron  faces  for  the  variables  with  three-dimensional 
dependency  may  be  viewed  as  a  polygon  side  in  a  two-dimensional 
space,  the  B ( I )  array  here  will  be  similar  to  that  in  Subsection  C.2. 

C.4  -  VARIABLES  WITH  HIGH-ORDER  SEARCH-RANGE  DEPENDENCY 

When  a  further  search-space  refinement  is  deemed  necessary,  such  that  the 
search  ranges  of  four  or  more  store  parameters  are  mutually  dependent,  there  is 
not  now  to  the  author's  knowledge  a  procedure  for  readily  defining  the 
constraint  hyperplanes.  Perhaps  the  best  ultimate  solution  would  be  a 
geometric  optimization  program  that  would  provide  the  arrangement  of  a 
specified  number  of  hyperplanes  that  minimizes  the  interior  volume  containing 
the  set  of  actual  store-inventory  points.  Until  such  a  capability  is 
available,  perhaps  the  best  approach  to  use,  when  it  is  important  to  consider 
lour  or  more  parameters  to  have  mutually  dependent  search  ranges,  is  to 
perform  searches  In  two  (or  possibly  more)  successive  three-parameter  spaces: 
The  three  parameters  that  are  judged  to  be  most  important  from  a 
flutter-speed  point  of  view  would  be  included  in  the  first  search  space,  and 
then  the  minimum-flutter-speed  point(s)  found  from  the  first  search  would  be 
used  as  starting  point (s)  for  a  search  in  a  second  space  that  would  include 
other  store  parameters  as  variables. 


C-5 


As  has  been  noted  in  Subsection  5.2.2,  one  special  case  of  a  1 our-pnrnmetor 
problem  that  can  be  readily  accommodated  at  present  is  the  situation  in  which 
pitch  and  vaw  inertia  are  considered  to  vary  in  an  identical  manner  (KCONST  - 
0  in  Item  8J )  .  The  number  of  variables  with  mutually  dependent  search  ranc.es 
is  now  reduced  by  one  from  a  program-logic  viewpoint,  so  th't  the  approach  of 
Subsection  C.3  is  directly  applicable.  Even  when  pitch  and  yaw  inertia  do 
vary  independently  to  a  significant  extent,  the  "slaving"  approach  will 
probably  constitute  a  good  first  approximation.  A  follow-up  search  can  then 
be  done,  ns  discussed  in  the  previous  paragraph,  in  a  space  containing  pitch 
and  vaw  inertias  as  separate  variables. 
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